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ABSTRACT 


A one dimensional, plug flow, non-adiabatic, heterogeneous steady state 
mathematical model is developed for tubular packed bed type exhaust catalytic 
converters for small spark ignition engine application. The proposed model incorporates a 
two step oxidation of propane and single step oxidation of carbon monoxide present in 
the exhaust stream with simultaneous diffusion and reactions on catalyst solid surface and 
in fluid phase. The mass and energy balance equations are coupled for solid phase and 
fluid phase reactions. The model predictions on conversion efficiency for carbon 
monoxide and hydrocarbon, temperatures and pressure drop are compared with a 
homogenous model and experimental results. 

Flow pattern in the muffler of a moped fitted with conical catalytic converter 
developed by M/s Indian Oil Corporation Limited, R&D Centre, Faridabad is studied 
using a CFD software. The catal5dic converter model is applied to this conical shaped 
converter. Predictions of model for the conical catalytic converter are compared with 
cylindrical catalytic converter of equal catalyst volume. 
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CHAPTER 1 


INTRODUCTION 


1.1 MOTIVATION: 

Many industrial activities, and particularly energy generation, industrial 
production and vehicular traffic have brought about a significant increase in 
environmental pollution. The rise in civilization is closely related to improvement in 
transportation of men and material. In the development of transport, the internal 
combustion engines, both petrol and diesel engines occupy a very important position. The 
petrol engine has provided rehable small power imits for personalized transport. The 
internal combustion engine is now known to be the major source of emissions of 
unbumed hydrocarbons, carbon monoxide and oxides of nitrogen into atmosphere.. 
These pollutants affect not only human health but also animal, plant and materials. 
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In the developing countries like India, importance of two-wheelers can not be 
neglected. The sale of new two-wheelers is approximately 10 times that of passenger cars 
and together they consume about 60% or more of the total gasoline in country [Pundir et. 
al. 1994]. The two-wheelers are highly polluting vehicles as most of these use two-stroke 
engine due to its simplicity, low cost and higher power to weight ratio compared to a 
four-stroke engine. The two-stroke engines produce a large amount of unbumed 
hydrocarbon and have poor efficiency due to its nature of operation in which short- 
circuiting of firesh mixture during scavenging and higher incidence of misfired 
combustion particularly at light load occur. The CO emissions are also high for the size 
of vehicle [O’ Sullivan, 1999]. But, NOx emissions are very low as a result of high 
residual gas fi'action, thereby lowering the combustion temperature [Hejwvood, 1989]. 
Two wheelers are the main contributors to urban air pollution in India particularly to CO 
and HC. In order to handle this problem, the Govt, of India has legislated for lower 
emission limits for all the mobile sources. 

In India emission levels were lowered in April 1996 and again in April 2000. 
Now, still tougher standards are being proposed for this decade [Pundir, 1997]. 

Research and development efforts have already started in order to meet the 
stricter norms in future. Various means have been suggested to control emissions for two- 
wheelers e.g. timing of intake and exhaust system, use of lean mixtures etc. But these 
have only given a marginal reduction in hydrocarbon and CO emissions. However, 
substantial emission reductions are possible only by elimination of scavenged through 
fuel loses or by after burning of unbumed hydrocarbon and CO in the exhaust system 
[Xiaobo et. al., 1996]. Since employing in-cylinder fuel injection is still on a laboratory 


2 



active catalyst surface area are sufficient to allow close to 100 percent conversion with 
high catalytic activity. 

1.2.1 Monolith Reactor System; 

This employs a ceramic or metallic monolith of honeycomb structure held in a 
metal can in the exhaust stream. The active catalyst material is impregnated into a hi ghl y 
porous alumina washcoat about 20 |am thick that is applied to the channel walls. The 
typical monolith has square-cross-section channel with inside dimension of ~ 1 mm 
separated by thin (0.15 to 0.3 mm) porous walls. The number of chaimels per square cm 
varies between about 30 and 60. The washcoat, 5 to 15 percent of the weight of the 
monolith, has a surface area of 100 to 200 m^/g. 

1.2.2 Packed Bed Reactor System: 

This system uses a bed of spherical ceramic pellets to provide a large surface area 
in contact with the flow. With pellet catalyst, the noble metal catalyst is impregnated onto 
the highly porous surface of the spherical alumina pellets to a depth of about 250 pm. 
The pellet material is chosen to have good crush and abrasion resistance after exposure to 
temperature of the order of 1000 °C. The gas flow is directed down through the bed to 
provide a large flow area and low pressure drop. 

The comparison of the packed beds and monolith structures shows, 

• In packed beds, attrition is a problem as a result of pellet motion during use. This 
problem may be reduced by using hard pellets. 
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• The use of low-density pellets to improve light-off characteristic may aggravate the 
problem of attrition. Also, as pellet loss occurs, voids are formed in the bed and 
b 5 ^assmg of exhaust gases may occur. 

• In packed beds, the gas flow is turbulent which results in high mass transfer rate 
while in monoUth structures, the gas flow is laminar which results in lower mass 
transfer rate. 

• Pressure drops are generally lower in monolith structures. 

In the present work, hard pellets i.e. of cordierite instead of alumina have been used. 

1.2.3 Oxidation Catalyst: 

The function of an oxidation catalyst is to oxidize CO and hydrocarbons to CO 2 
and water in exhaust gas stream. About half the hydrocarbons emitted by the SI engine 
are unbumed fuel components. The saturated hydrocarbons (which comprise 20 to 30 
percent) are most difficult to oxidize. The ease of oxidation increases with increasing 
molecular weight. Sufficient oxygen must be present to oxidize the CO and HC. This 
may be supplied by the engine itself running lean of stoichiometric or by injection of 
secondary air upstream of the converter. Becaixse of their high intrinsic activity, noble 
metals are most suitable as the catalytic material. They show higher specific activity for 
HC oxidation, are more thermally resistant to loss of low-temperature activity, and are 
much less deactivated by the sulphur in the fiiel than base metal oxides. A mixture of 
platinum (Pt) and palladium (Pd) is most commonly used. Pure noble metals sinter 
rapidly in the 500 to 900 °C temperature range experienced by exhaust catalyst. 
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Indian Oil Corporation Limited with The Associated Cement Cos. Ltd. (ACC) has 
developed an alternate material for pellet beads, which has good properties over 
conventional alumina pellet beads. 


1.2.4 NO Catalyst: 

NO is removed by reduction reactions using the CO, hydrocarbons, and H2 in the 
exhaust. Alumina-supported noble metal catalysts reduce NO with CO-H2 mixtures. The 
NO-reduction activity of noble metals is in the order Ru>Rh>Pd>Pt. Ruthenium (Ru) 
and Rhodium (Rh) produce considerably less NH3 than Pd or Pt under slightly rich 
conditions. While these properties make Ruthenium a desirable NO catalyst, it forms 
volatile oxides under oxidizing conditions which result in loss of Ruthenium from 
alumina support. Therefore Rh is generally used for reduction of NO. 

1.2.5 Three-wav Catalyst; 

If an engine is operated at all times with an air/fuel ratio at or close to 
stoichiometric, then both NO reduction and CO and HC oxidation can be done in a single 
catalyst bed. Enough reducing gases CO and H2 will be present to reduce NO and enough 
O2 to oxidize the CO and hydrocarbons. Such a catalyst is called ‘Three-way Catalyst’ 
since it removes all three pollutants simultaneously. There is a narrow range of air/fuel 
ratio near stoichiometric in which high conversion efficiency for all three pollutants are 
achieved. The width of this window is narrow, about 0.1 air/fuel ratio for catalyst with 
high mileage use, and depends on catalyst formulation and engine operating conditions. 
Because of cyclic variations in exhaust gas composition about a set point close to 
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stoichiometric, it is desirable that catalyst be able to reduce NO when a slight excess of 
oxygen is present (on the lean side) and remove CO and HC when there is a slight 
deficiency of oxygen (on the rich side). To remove NO under slightly lean-of- 
stoichiometric condition, the catalyst must react the CO, H2, or HC with NO rather than 
O2, as the exhaust gas passes through the catalyst bed. Rhodium is the principle 
ingredient used in commercial catalysts to remove NO. It is very active for NO reduction, 
is much inhibited by CO and sulfur compounds, and produces less NH3 than Pt. 

1.2.6 Catalytic Converter Design Considerations: 

The requirements of catalytic systems, in addition of the main task of achieving 
legislative emission requirements, are: 

• Minimal space requirement, 

• Minimal additional weight, 

• Limited complexity, 

• Minimal pressure drop, 

• Minimal additional cost , 

• Suitable production technology for mass production and 

• High durability. 

Catalytic oxidation reactions being highly exothermic in nature, catalytic reaction 
steps giving rise to temperature gradients. The released energy is first transported by 
conduction within the catalyst and then by convective heat transfer to the bulk fluid. If 
sufficiently high temperatures are encountered within the reactor, homogeneous gas 
phase reactions also contribute to the overall energy balance. As the gases flow past the 
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converter, depending on its geometry, substantial pressure drop may be encountered 
which may severely affect the engine performance. The catalyst support system should 
permit successful completion of the above processes and achieve the desired results. 

1.3 MODELING OF CATALYTIC CONVERTERS: 

A mathematical model is a system of equations, which describes the relationship 
among the various physical and chemical variables governing the process. 

Mathematical models are classified into several types. One extreme is the empirical 
model, in which experimental data are described by an arbitrary set of fimctions, the only 
criterion being the ‘goodness of fit’ determined by statistical means. An experimental 
model may describe a set of experimental data and be able to predict new values of a 
dependent variable within the range of experimentally measured parameters. The danger 
in using such a model is tendency to extrapolate to regions outside of the original 
ejqperimental parameters. On the other hand, a purely empirical model should not 
necessarily be discarded simply because it has little or no mechanistic base. 

The other type is phenomenological model, which is based on the fundamental 
physio-chemical process encoimtered. If it is possible to derive a theoretical model, 
which describes all such processes, then the model can be used directly for design 
without the need for experiments. This is invariably an impossible task for real systems, 
usually because some physical constants must be measured experimentally. The 
phenomenological model can be divided into ‘Lumped Parameter Model’ and 
‘Distributed Parameter Model’. A lumped parameter model ignores some or all of the 
spatial variations, which may exist for one or more of the properties of the system and 
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uses an average value. A distributed parameter model accounts for all spatial property 
variations. 

A practical and useful approach may be to derive -as scientifically meaningful a 
model as possible with a few constants as necessary, and then to determine each constant 
in a separate experiment. The final model should be tested for its ability to predict the 
behaviour of the full process. 

Factors, which affect the selection of the type of model, include the following: 

• The precision to which the model parameters are known, 

• The desired accuracy and precision of the solution, 

• The importance of various phenomena, 

• The power of the computer. 

Enhanced computational capabilities have changed the method of design and 
development so contribution of mathematical modeling gained importance while earlier 
works were based on experiments. 

There are four distinct areas for the modeling of catalytic converter: 

• Chemical Reaction Kinetics, 

• Heat Transfer & Thermodynamics, 

• Gas Dynamics, 

• Diffusion and Adsorption. 

Though the theories have been developed but their successful applications as an 
usable engine design tool is rather sporadic in nature. While mathematical modeling of a 
complete engine is now a standard tool for the engine designer, modeling of the catalytic 
converter systems has not evolved at a comparative pace. 
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1.4 LITERATURE SURVEY: 


A number of studies on mathematical modeling of the catal 3 dic converter have 
been carried out since early 1970’s both on packed bed and monolith catalytic converters. 

1.4.1 Packed Bed Catalytic Converters: 

Kuo et. al. [1971], four years before the implementation of the first emissions 
reduction legislation in the USA requiring the use of catal 5 d:ic converters, developed 
mathematical model that predicts the performance of cataljhic converter system. The 
model contains all the important system parameters and process variables e.g. catalyst 
density and heat capacity, catalyst activity, initial bed temperature, inlet exhaust gas 
composition and converter location. 

Se et. al. [1980] developed a detailed mathematical model of a single catalyst 
pellet to describe its behavior under transient conditions encountered during the warm-up 
period of automobile exhaust catalytic converters. This model is capable of describing 
the time behavior of a composite catalyst pellet, and thus provides a convenient means of 
examining the effects of poisoning penetration and various noble metal impregnation 
strategies. They found that catalyst light-off was delayed with increasing carbon 
monoxide or propylene concentration in the exhaust gas, whereas increasing the oxygen 
concentration improves both light-off and steady state activity and faster light-off occurs 
when a catalyst pellet exposed to a higher exhaust gas flow rate. 

Oh et. al. [1985] developed a transient mathematical model which described the 
dynamic behaviour of packed bed catalytic converter during warm up. They found that 
the model prediction agreed very well with the results of engine dynamometer 
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experiments for three Pt-alumina catalyst of different properties and can be used to 
predict the warm-up performance of packed bed catalytic converter from specified 
system parameters and inlet exhaiist gas conditions. Oh et. al.[1985] further found that 
the model could also be applied to dual bed converter because the latter converter 
normally operates in the oxidizing mode during warm-up. 

Khandekar S. [2000] developed a steady state mathematical model of packed bed 
catalytic converters for two-stroke petrol engines. He used homogeneous model in which 
it was assumed that concentrations and temperatures of species at pellet surfece and in 
gas phase are same. The model predictions were compared with experimental results. A 
new catalytic rate of reaction for carbon monoxide and propane was proposed which is 
being described in section 1.4.3. He further concluded that a heterogeneous model is 
needed for improvement in prediction ability of the model. 

1.4.2 MonoKth Catalytic Converters: 

Larry et. al. [1976] applied orthogonal collection method to solve the 
mathematical model of monolithic catal5dic converter. The orthogonal collection method 
was developed for application to problems having irregular geometeries. Most of the 
previous works were on rectangular region. They also included transient heat and mass 
transfer in three dimensions. They found that more accurate results were obtained for the 
asymptotic Nusselt number for the geometeries like square, rectangles with aspect ratio 
of 0.5 and 0.25. 

Heck et. al. [1976] developed mathematical model which account for 
simultaneous heat transfer, mass transfer and chemical reaction in the oxidation of carbon 


11 




monoxide over platinum containing monoliths. A two-dimensional model was used to 
predict unusual behavior of the Nusselt number in the presence of rapid reaction. It was 
found that monolithic reactor behavior was dominated by the position of the light-off 
point, where the catalytic walk temperature took an upward leap, the wall concentration 
fallen to zero, and the reaction becomes mass transport limited. 

Naoki et. al. [1996] analyzed numerically the conversion characteristics of 
monolith exhaust catalyst under warming-up conditions and verified it with experimental 
results. They found that there is large maldistribution of species concentration inside the 
catalytic converter during warm-up. In the outer region of monolith, CO and HC pass 
through the channel without being converted. They also found that during warm-up, 
about 50% of the total HC mass in catalyst out-gas are exhausted firom the channels in the 
outer 25% of the cross-sectional area of the converter. Conversion performance of the 
exhaust monolith during warm-up is greatly improved by both heat insulation and high 
noble metal loading on the outer region of the upstream part of the monolith. It was seen 
that decreasing metal loading on the downstream part of the monolith did not reduce the 
conversion and high noble metal loading on the inner region of the upstream part of the 
monolith accelerated the conversion rates at an early stage of warming-up process. 

1.4.3 Chemical Kinetics: 

For modeling of the catalytic converters, reaction rate constants for oxidation of 
carbon monoxide and unbumed hydrocarbons are required. 

Sterling et. al. [1973] studied kinetics of carbon monoxide and propylene 
oxidation on a platinum- alumina catalyst with synthetic gas mixture between 205 °C and 
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372 °C. The rates of oxidation of carbon monoxide and propylene increase with 

increasing concentration and are inhibited by carbon monoxide, propylene and nitric 
oxide. They formulated complex kinetic equations are of the form, 


KXcdl02) 

R[e) 

r = 

where kri and kri are rate constants and R(9) is a resistance term which includes 
the inhibition effect of carbon monoxide , propylene and NO on reaction rate. 

Hegedus et. al. [1977] used an integral reactor to study the CO-O 2 reaction over 
Pt-coated beads. They presented a rate expression in the same form as Voltz et. al. [1973] 
without the inhibition terms for propene (CsHe) and NO. They obtained a different 
activation energy (about 40% lower) and CO adsorption equilibrium constant (5% 
higher) than Voltz et. al. [1973] and attributed the difference to the influences of CsHe 
and NO. Also, stable isothermal multiplicity was observed which was investigated in the 
conversion-temperature-mass flow domain. 

Cant et. al. [1978] studied steady state oxidation of carbon monoxide over 
supported noble metals with particular reference to platinum. Activation energy, kinetic 
orders and relative activities had been determined for the oxidation of CO by O 2 over five 
supported noble metals. The results were compared with the prediction of a model, which 
used parameters derived from measurements made with only reactant present in the gas 
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phase. The results showed much lower activation energy and only a shght negative order 
dependence on CO. 

Herkowitz et. al. [1983] studied CO oxidation on platinum supported catalysts. 
The conversion of carbon monoxide to carbon dioxide was measured in a recycle reactor 
automatically controlled by a computer. Two platinum/Si02 catalyst of different 
dispersion were employed. The rate data were correlated by three LHHW kinetic models. 
Which are listed in table 1.1. 


Table 1.1 


Model 

Mechanism 


1 

CO + soCO, 
O2+25O 20, 
CO, + O, — > CO 2 + 2s 


II 111 

CO+50CO, CO + sc^COs 

O^-^soO^s O^ + soO^s 

COs + O2S CO2 +• O, + 5 COs + O2S C 0 ^S 2 

(2Q2S2 + CO, — 2CO2 + 3 s 


Kinetic Expression, r 

K ' ^co ' 

j^l/2 ^1/2 

^02 ' ^CO * '-'Os 

K ‘ ^CO ' ^Oz ’ ^ CO ' ^Oz K * ^CO 

’ ^Oz * ^CO ‘ ^Oz 

{i + K^o-Cco+K'<'J 

i^ + Kco-Cco+Ko^-Coy {l + K,o- 

^ CO ^Oz * ^Oz f 

Paramters 

1 

2 

1 2 1 

2 

A^, Kg mol I{kg.surf,Pt) 2.8 x 10* 

9.1x10* 

8.1x10^ 4.8x10^ 1.3x10' 

7.9x10^ 

A 2 ,m^ 1 Kg mol 

39 

37 

16 13 23 

17 

A^,m^ I Kgmol 

1.8x10'^ 

3.0x10"^ 

2.4x10* 1.3x10* 2.0x10* 

1.1x10* 


7680 

8000 

7980 


0^2, K 

-2210 

-3230 

-3330 


ay,K 

1120 

-230 

-310 



a, s indicate an active site, 

b. k, = A, ■ expf- = 4 - ejqpf- Ko^ = 4 ■ expj 
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Khandekar S. [2000] used LHHW model with modified catalytic rate of reactions 
for carbon monoxide and also propane. For catalytic rate of reaction for carbon 
monoxide, he modified activation energy term as proposed by Hekowitz et. al. [1983] 
from 7980 to 8300. Thus, 
k'-k -k 


■Rco = 


where 


'CO 


(l + kf~Q ■ C(2Q + • Cq^ ) 


k' = A^- exp 


-8300 


kco=A^ - exp 


V T j 
^- 3330 ^ 


^ 3 ■ exp 


k T 
•310' 


/ 


where A^,A2 and A^ are constants. 

For propane oxidation the pre-exponential constants in the catalytic rate of 
reaction as proposed by Wanke[1973] was modified from 3.15*10"^ to 12.8*10'^. 


Thus 


=12.8x10" - exp 


^ 89126^ 


V 




where Rg is universal gas constant. 


A brief summary of the various forms of the rate of reaction for carbon monoxide 
is given below in Table 2.2: 
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Table 2.2 


Rate Law Form 

Proposed By 

Reference 

nC _ Oi 

^CO ■” yO.5 
^CO 

Shishu (1972) 

(Boehman etal. 1 993) 

Jr •¥ Y 

pC __ 

(Yo^+kJ^o) 

Young (1974) 


k-YaJco 

^CO - . X 2 

Shishu (1972) 


k-YoYco 

Hegedus et.al. (1977) 


k- KqqKq^ 

^CO / \2 

{i + K,oCco+KoCoJ 

Herkowitz et.al. (1983) 



Montreuil et.al. (1992) 



1.5 STATREMENT OF THE PROBLEM: 

For the purpose of present study, the overall domain of the problem has been 
subdivided into the following parts: 

I. To develop a one-dimensional, plug flow, non-adiabatic, heterogeneous, steady 
state mathematical model for tubular packed bed exhaust catalytic .The model should 
predict the steady state conversion, temperature and pressure profiles along the axial 
direction. 
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III. The model should incorporate simultaneous diffusion and reaction phenomenon 
occurring in the catalyst washcoated cordorite pellets. 

IV. To compare the model with homogenous model and experimental values. 

V. To validate the model with experimental data. 

VI. To study flow field in the exhaust muffler of a converter designed for a moped 

using CFD software. 

VII. To carry out parametric analysis with the mathematical model for above moped 
catal3d;ic converter. 
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CHAPTER 2 


MATHEMATICAL MODEL 


2.1 INTRODUCTION: 

The heterogeneous mathematical modeling of the packed bed catalytic converter, 
involves the formulation of mathematical equations describing the following: 

(al The Basic Governing Equations: 

• The Solid Phase Mole Balance Equation. 

• The Fluid Phase Mole Balance Equation. 

• The Momentum Balance Equation. 

• The Solid Phase Energy Balance Equation. 

• The Fluid Phase Energy Balance Equation. 

(bl The Chemical Kinetics; 

• The Rate of Homogeneous Reaction for Carbon Monoxide. 

• The Rate of Homogeneous Reaction for hydrocarbons, in this case. Propane. 

• The Rate of Catalytic Reaction for Carbon Monoxide. 
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• The Rate of Catalytic Reaction for Propane. 

(cl The Physical Properties: 

• Gas Viscosity Data. 

• Specific Heats for all the species involved. 

• Heat of Reaction for Propane and Carbon Monoxide oxidation. 

• Effective Diffusivity of the Catalyst Pellet. 

(dl The Heat and Mass Transfer: 

• The Thermal Conductivity of the Pellet. 

• The Thermal Conductivity of the Fluid. 

• Solid Surface to Fluid Heat Transfer Coefficient. 

• Bed to Wall Heat Transfer Coefficient. 

• Overall Heat Transfer Coefficient for the Packed Bed. 

• Fluid to Solid Surfece Mass Transfer Coefficient for all the species involved. 

(el The Pellet Diffusion and Effectiveness Factor; 

• The Catalyst Pellet Effectiveness Factor for Carbon Monoxide. 

• The Catalyst Pellet Effectiveness Factor for Propane. 

2.2 THE BASIC GOVERNING EQUATION: 

In this section, the fundamental governing equations in the present study are 
derived. Control volume as shown in Fig. 2.1 is chosen for which mole, momentum and 
energy balance equations are applied for solid and fluid phase. Fig. 2.2 shows the process 
of mass transfer from bulk flow to the solid surface of the pellet and heat transfer from 
solid surface of the pellet to the bulk flow. 
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2.2.1 The Solid Phase Mole Balance Equation: 


The rate of disappearance of species in the solid catalyst is; 


Moles 


Moles 

transported to 


reacted in 

catalyst surface 


catalytic reaction 


which can be expressed in terms of the incremental surface area as: 


k,yAS-c,-{Y' -Y;)=n-(-R;)-Av-p,ii-^) ( 2 . 1 ) 

where is the fluid to solid mass transfer coefficient for the ‘z’th species, Cf is the 
molar density (concentration), Y( and Yf are mole fractions of ‘i 'th component at fluid 
and solid phase respectively, t] is the effectiveness factor, (Ri)c is the catalytic rate of 
formation of V ’th species, pc is the mass density of catalyst pallet, <j) is the porosity of 
the catalyst bed. 


Noting that. 


a 


V 


^V 


( 2 . 2 ) 


The equation can be simplified to: 

■ cr (y - j;' )- n ■ (- )■ Pc ■ (1 - <>) = 0. 


(2.3) 


where, Uv is the solid particle surface area per unit bed volume. 
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2.2.2 The Fluid Phase Mole Balance Equation: 


The mole balance for the fluid phase is, 


\Moles in] - \Moles out\ ■ 


Moles reacted in 
homogeneous 
reaction 


Moles transported 
to catalyst surface 


= 0. 


Which is, 

\ 




(2.4) 


Noting that F,- = Vs-Ac- Q = Vs ■ Ac ■ Cf Yi and dividing by AV =Ac-Az , introducing 
and as the volume element, final mole balance equation is. 


az 


(2.5) 


where F,- is the molar flow rate, is the superficial velocity and i? * is homogenous rate 
of formation of species ‘i ’ ,Acis the cross sectional area of packed bed. 

(a) Multiple Reactions: 

In the present study, multiple reactions are modeled. A two step oxidation 
reaction of propane, which is first oxidized, to carbon monoxide and then the carbon 
monoxide later, gets converted to carbon dioxide. Detailed explanation of this two-step 
model is presented in section 2.3. 

Necessary modifications, as described below, have to be incorporated into the 
governing equations for accounting for the multiple reactions. The general hydrocarbon 
reaction for two step oxidation is given by: 
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(2.6-a) 

>^C0{g) +—02{g) -^nC02{g) 

Step II 

(2.6-b) 

Taking «=3 and m=8, we get, 



C2H,{g) + 3.502{g) -» 3C0{g)+4H20{g) 

Step I 

(2.6-c) 

cofe)+io,(«) -^CO,{g) 

Step II 

(2.6-d) 

Since we explicitly know the intrinsic rates of reaction for propane (CsHg) and 

carbon monoxide (CO), we can apply the stoichiometric relationship between the reacting 

species to get the following relationship for the rate 

of reactions for both catalytic 

reaction and homogeneous reaction, from Equation 2.6-c, 





(2.7-a) 

(2.7-b) 




• ^("”^C3//g ) 


(2.7-c) 



(2.7-d) 

and from Equation 2.6-d, we get: 



-rco =l-0(-^co) 


(2.8-a) 

-'o, =0-5(-^co) 


(2.8-b) 

~^C02 • 


(2.8-c) 
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The subscripts / and II refers to the reaction steps. Thus applying the solid phase 
mole balance Equation 2. 3, we get, 



(2.9-a) 

c, -t. -o.-fe -i'i)=--T{3.5(''c:«.)+0.5{''S)}'Pc'(l-'^) 

(2.9-b) 

'^v '{^co -(l — 

(2.9-c) 

CfK-^.- fe - ) = -1 ■ {- ‘‘.ofe. )}■ Pc ■ (1 - ■(■) 

(2.9-d) 

Cc ■ i. • a. ■ fe. - %. ) = -1 • !- 1 -ote ))■ Pc ■ (1 - ) 

(2.9-e) 

and applying the fluid phase mole balance equation 2.5, we get. 

In general, for catalytic combustion, the molar flow rate will be changing as the 

reaction progresses. As only equal number of moles occupy equal volumes in the gas 

phase, the volumetric flow rate will also change. Thus, we get. 


(2.10-a) 

— )+ 0.5te)}-^. -k.-a,-C, ( y^ - ) 

(2.10-b) 

= {- s-ofe. )+ 1 .ofe ))• - 1. ■ a, ■ (r/c - Ki ) 

(2.10-c) 

= -4.ofe. ) ■ ^ fe -yi,o) 

(2.10-d) 

-a, ■c,(Yy„^ -rij 
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(2.10-e) 



(2.11 -a) 


^TOTAL — -^Cj/Zg ^CO ^CO-^ 

^TOTAL ~ ■*■ -^O, “*■ ^CO “*" ^Hi,0 

where, 


+ F„ 


+f:. 


(2.1 l-b) 


Ftotal ~ Total molar flow rate at a point, and = Entering total molar flow rate of 

the species. 

The next step is to replace all the concentrations terms in the rate laws in terms of 
molar flow-rates. For a flow system, the concentration of species at a given point can be 
determined from the molar flow rate and the volumetric flow rate of the specie. We know 
that. 


r P 

V^ijnr 


(2.12-a) 


e R,-T, 


F P 

n — ^ Tm _ 


Qo R,-T^n 


(2.12-b) 


where, 

Ct = Total concentration, P= Absolute pressure, i?g=Universal Gas Constant, r=Absolute 
temperature and Q= Volumetric Flow Rate of species, Tf is fluid phase temperature. The 
subscript in refers to the inlet conditions. 

From Equation 2.12-a, b, we get. 


Q = Qi. 




\^rm ) 








z. 


(2.13) 


Thus, in general. 
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(2.14) 


where Cj refers to the concentration of species j . 


2.2.3 The Momentum Balance Eauation: 


For Newtonian fluids, the one-dimensional form of the Darcy’s law for flow 


through porous media, giving the linear relationship between velocity and pressure drop 
is given by [Bear, 1972]: 


dz K 



= ^ • V 


(2.15) 

(2.16) 


The above law is valid at low flow rates. At higher flow rates fluid inertial effects 
become important. Therefore, an extension to Darcy’s law, known as the Forchheimer 
equation [Hayes and Kolaczkowski, 1997] was proposed to account for those effects. The 
general form is given by. 


I C-p-Vs 
dz Ik 4k 


(2.17) 


The constant C may assume various values as per the actual working situation. A 
commonly used permeability for packed beds of uniformly sized smooth particles is 
obtained from the Black-Kozeny friction factor analysis, supplemented with experimental 
results [Bird et. al., 1960], which gives. 


Dl ■<!>' 


(2.18) 


The equivalent particle diameter, is defined as. 
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D,=- 


a,. = 


where, 

«v 

particle surface area to volume 
particle volume 


(2.19) 


For the case of spherical pellet = actual diameter of the pellet. A commonly expressed 


form of the constant C is given by, 
0.143 


C = - 




1.5 


( 2 . 20 ) 


Substituting Equations 2.18 and 2.20 in 2.17 we get, 
l50-u-(\-(j)) 


dP v.-(l-<^) 

dz 




D. 


( 2 . 21 ) 


The above equation is called the Ergun Equation [Ergun, 1952] widely used for 
predicting pressure drops in packed beds. The value of the constants appearing in 
Equation 2.21, viz. 150 and 1.75 was originally proposed. The values depend on the 
particle surface roughness, particle size distribution, particle shape and packing 
homogeneity. 

In catalytic combustion reaction the fluid is a compressible gas. Thus the fluid 
density depends on the temperature and the pressure. The equation of continuity can be 
applied to one-dimensional flow in a packed tube of constant cross-sectional area to get: 
d{p-Vs), 


dz 


• = 0 or p • v, = p.„v^^ = G = cons tant 


( 2 . 22 ) 


which modifies Equation 2.21 to. 


dP _ G-{l-(l>) 


dz 


pD-^- 


150m(1-<>) 


D. 


+ 


(1.75 -G) 


(2.23) 


Equation 2.23 can be further modified by using. 
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p 


(2.24-a) 


Cr = 


Ft 


Q RT,z 


where, Z and Z,.„ are the compressibility factors, 


Fr- 

C _ Tin 

Ti'm 


"Tin 


RTnZin 


Q = Q. 


and 


P« =P 


f FtYkY rA 


K^ToA ^ J 


\Rinj 


(neglecting the ratio of compressibility factor ). 




J 


( x> \ 


f T \ 


J 


\RinJ 


to obtain, 
dP _ G •(!-(/>) 


Pin- Dp Y 




f TP p \ 


J 


(tA 


yZinj 


(2.24-b) 

(2.25) 


(2.26) 


(2.27) 


This is the final form of the momentum equation for the plug flow tubular reactor. 

2.2.4 The Solid Phase Energy Balance Equation: 

The solid phase consists of the catalysts, which is not flowing, and the gas, which 
fills the pore volume of the catalyst. We can consider the solid as an open system with 
diffusion of mass in and out. The gas reacts and heat is released. At steady state this must 
be balanced by heat that is transferred from the solid. In other words, 

[Heat released by reaction] = \Heat transfer from solid] 
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In this analysis heat is transferred from the solid by convection. Heat transfer by 
conduction and radiation are neglected and the whole pellet is assumed as having the 
uniform temperature. 


The heat released by reaction is. 


\Heat releaed by reaction\ = -(l - (^) • p,- • 77 • (- i?/ )• AH ^ 

where AHr is the heat of reaction at solid phase temperature and hfs is the heat transfer 
coefficient between fluid and catalyst surface. 


The solid phase energy balance equation is therefore, 

(1 -'(')■ Pc ■ t) AHJ + V, ■ O. ■ (r, - )- 0 

2.2.5 The Fluid Phase Energy Balance Equation: 

The steady state fluid phase energy balance equation is. 


(2.28) 


Enthalpy increase 


Heat added 


Heat generated 


Heat added 

as a result of 

= 

togas 

+ 

in gas phase 

+ 

from 

temperature rise 


from surface 


reaction 


surrounding 


The left-hand side is. 




(2.29) 


where Qv is volumetric flow rate, pis mass density, Cp is constant pressure heat capacity. 
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The area of heat transfer in the elemental volume is the circumference of the channel 
multiplied by the length, that is, 

bAj=K-Dj-ba ( 2 . 30 ) 


where Dt is the reactor diameter. The heat exchange with surroundings depends on an 
overall heat transfer coefficient, which includes the heat transfer at the inside reactor 
wall, the outside reactor surface and conduction through the wall, and the driving force, 
which is the temperature difference between the fluid in the reactor and the surroundings. 
The heat loss from the elemental volume element is then. 


Heat added 
from surroundings 


= q 


conv 


:=n-D,-Az-u\T,^^-Tf) 


( 2 . 31 ) 


where is the temperature of the heat transfer fluid surrounding the reactor. 
The right hand side is. 




( 2 . 32 ) 


where AH/ is heat of reaction at fluid phase temperature and hfs is the pellet to fluid heat 
transfer coefficient. 

Therefore energy balance over AV is, 

-Q^-p-C,-AT^-AV-<l>\-R^\lsHi+AS-h^,\T,-Tf)+7t-D^-Az-U-[T,^^-T^)=0 

(233) 
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Dividing through hy AF = Ac ■ Az, noting again that, 

AS O K t 

= , V, and AF = --Dj. -Az (2.34) 

' AV ' Ac A ^ ^ 

gives, 

-V, -p-C, •(-«*)♦+*» -a. (2-35) 

Az Dj 

Taking the limit as AV~*' dv, we obtain the ordinary differential equation, 

dT A. 

-v,-p-C,-^-AH.-(-i!,*).^. + ft„-a.-{r,-r^)+— (2.36) 

The temperature at the reactor inlet must be known to solve this IVP. 

2.3 THE CHEMICAL KINETICS: 

The rate expression, essentially a mathematical model of the reaction, relates the 
rate of reaction to the local temperature, concentration of the reactants and two or more 
empirical constants. When expressed in terms of a local value of temperature and 
concentration, it is termed as the intrinsic rate expression. The global rate is a rate that is 
observed in a reactor, and may include much heat and mass transfer effects combined 
together. The reaction mechanism is a description of all the steps, which occur before the 
reactants are turned into final products. The development of reaction rate expressions is 
often quite difficult, drawing upon empirical evidence as well as theory. 

The rate expressions for elementary reactions of the form, oA+bB — ^^cC+ dD are 
almost always expressed in the Arrhenius type model, 
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where ki is called the pre-exponential constant, Ea = activation energy, Rg is the 
universal gas constant and a and /3 are real numbers. It is difficult, if not impossible, to 
predict the pre-exponential factor and the activation energy theoretically and these values 
are usually determined from experimental data. 

$ 2.3.1 Homoeeneous Oxidation of Carbon Monoxide: 

Apart from its presence in the exhaust gas, carbon monoxide appears as an 
intermediate chemical species in the combustion of hydrocarbons. When hydrocarbons 
are oxidized, CO appears first and is later converted to CO 2 in a slow secondary reaction. 
The overall reaction may be represented as: 


Stsp I HC O 2 — + 
Step II : CO + ~ O2 — ^ CO 2 


(2.38-a) 


(2.38-b) 


While there have been many studies of the oxidation of CO, a frequently used rate 
expression for homogeneous CO combustion is [Hayes et. al. 1997], 


=1.26e’°exp 


^ 20131^ 


y J 


^ ^ 0.25 ^ 0.5 

^CO 




(2.39) 


where concentrations are in moL/m^ and temperature is in Kelvin. 
2.3.2 Homogeneous Oxidation of Propane: 


Much of the work done in determining rate equations for the homogeneous 
combustion of hydrocarbons has been based on the study of laminar flame propagation. 
The actual mechanism of hydrocarbon combustion is very complex and the inclusion of 



this mechanism in a reactor model will yield a very expensive computational problem. 
Furthermore, for most catalytic combustion applications it is probably not necessary to go 
into such details. 

The simplest model has a single step. This model contains only four species and no 
intermediate concentration is included. 

(g) + (f + (g) n CO(g) + fH20ig) (2.40) 

Some success at predicting flame speeds has been shown using one step models, 
but they obviously have some drawbacks. One of these is that they predict the immediate 
formation of CO 2 in one step, whereas it is known that CO is produced first, and then 
reacts relatively slowly to form CO 2 . To account for these shortcomings various multi- 
steps models have been proposed. However it must be realized that for every new species 
introduced in the homogeneous reaction scheme, a corresponding catalytic rate 
expression must be available. Keeping this in mind it is probably best to limit the model 
choice to at the most a two step model for the homogeneous oxidation of propane. 

Therefore for the present study, a two step model introduced by Dryer, Hautman 
et. al. [1981] and Westbrook [1984] has been incorporated. The overall reaction is 
summarized as, 

C 3 H 3 + 3 . 5 O 2 -> SCO + AH^O (2.41 -a) 


C0 + -0, ^CO, 
2 ' ' 


(2.41-b) 
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The CO oxidation step is same as given in the previous section. The rate 
expression for homogeneous combustion of propane gas as given by [Hayes, 1 997] used 
in present study is. 



where all concentrations are in mol/m^ , temperature is in Kelvin and the rate is expressed 
in mol/m -s. 


2.3.3 Catalytic Oxidation of Carbon Monoxide: 


The Langmuir Hinshelwood Hougen Watson (LHHW) reaction model is 
frequently used to model the catalytic oxidation of carbon monoxide. This model 
assumes that the rate determining step (i.e. that step which is intrinsically slow) in the 
overall mechanism of CO oxidation is the reaction between adsorbed CO and adsorbed 
oxygen atoms or molecules [Scott Fogler, 1994]. The common feature of all LHHW 
models is that the rate equation is expressed as a ratio. The numerator is referred to as the 
kinetic driving force. The denominator represents the inhibition force due to adsorption 
of the species on the surface. The general form of LHHW models for catalytic oxidation 
of CO is represented as, 


k' 

(l + ^CO^CO + ^02 ^ 


(2.43) 


It can be seen that, at any given temperature, an increase in reactant concentration 
will result jp an increase in the kinetic driving force, which has the tendency to increase 
die rate of reaction. However, this will also increase the value of the adsorption inhibition 




term that will tend to decrease the rate of reaction. Thus, the net result solely depends on 
the local conditions. 

LHHW rate equations are also well known to give rise to multiple steady states 
[Nauman, 1987]. From mathematical point of view this means that for any given set of 
inlet conditions to the reactor, there exists more than one set of solution variables that 
satisfies the set of governing equations. The steady state at which the reactor will actually 
function depends on how operating conditions were varied, or in other words the history 
of the reactor. Plug flow reactors lack any internal mechanism for memory. There is no 
axial dispersion of heat or mass. What has happened previously has no effect on what is 
happening now. This fact prevents steady state multiplicity in plug flow reactors, unless 
there is some form of external feedback. Given a set of operating conditions, only one 
output is possible. 

Various forms of the above model, as reported in the literature have been 
incorporated and tested. 

2.3.4 Catalytic Oxidation of Propane: 

Propane is a typical compound frequently used as a representative to model the 
fast burning components of an automobile exhaust, as reaction rate constants for typical 
gasoline are not available. 

Typical rate expressions encountered in the literature are presented in the form of 
power law models. Wanke [1973] proposed the intrinsic rate model used in the present 
study. The study was based on oxidation of propane in air at atmospheric pressure over a 
commercial catalyst. The spherical beads have the active catalyst deposited in the 
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washcoat near the surface of the pellets. The rate expression was reported to be linear 
having the form, 


-R^ch =3.15e"exp 


89126 


V J 


C. 




(2.44-a) 


The rate has the units of mol/kg-s, the mass being the total mass of the entire pellets 
including the weight of the substrate. The above model does not take into account the 
inhibition effects, if any, by other chemical species present. 

Khandekar [2000] modified it with new value of pre-exponential coefficient, 
which has been used in the present problem. This is following. 


-RIj, =12.85x 10" exp 


^ 89126^ 


V j 


C. 


C.H, 


(2.44-b) 


2.4 PHYSICAL PROPERTIES: 


The physical property models, required for the overall model development are as 
follows: 

2.4.1 Gas Viscosity: 

The viscosity of the gases is approximated to the viscosity data for air given in 
Holman [1997] and is modeled as. 


/I = 3.1227e-^ +5.47613e-*r-1.93284e-"r'-t-3.65582e-''r" (2.45) 


where T is in Kelvin and }jl is in Pa.s. 

2.4.2 Specific Heat of Species: 

The specific heat of the species, involved in the chemical reaction, as a function 
of temperature are modeled as [Hayes, 1997] 
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Cp {C,H , ) = -4.04 + 30.43e" r - 1 5.70e-'r' + 3 1 
C^iCO^) =22.22 + 5.971 e-'r-3.495e''r' +7.457e'Y' 
C^{CO) = 28.11 + 0.067e‘'r+0.536e-'r' -2.218e'Y' 
C/O 2 ) = 25.44 + 1 .5 1 8 e" Y - 0.7 1 5e' Y' +1.31 0e‘ Y' 

C/TVj) = 28.85 -0.157e-Y +0.807e'Y' -2.868e'Y' 

where T is in Kelvin to give Cp in Joule/mol.K. 

2.4.3 Heat of Reaction; 

For the two step reaction represented eis, 

Cji/g + 3 . 5 O 2 -^3C0 + 4H,0 
CO + 0.5 O 2 — ^ CO 2 

by definition, we know that, 

jAC,-<irand AC, = ^ j.C, - 

298.15 products reac tan ts 


(2.46-a) 

(2.46-b) 

(2.46-c) 

(2.46-d) 

(2.46-d) 


(2.47-a) 

(2.47-b) 


(2.48) 


substituting the values from Equation 2.44, we get, 

(aC„ ) = 128.09 - 0.344734r + 24.0253e“Y^ + 56.2985e’Y^ 

(aC^)^^ = -18.61 + 0.050457’ -3.6741e"Y^ +9.02e-Y" 

and Equation 2.45-a transforms as, 

A/7^.c,//,(^) = -1220779.-128.ir + 0.172r^ -S.OOSe’Y' +14.07e-Y'* 
^r,co{t) = -279380. -18.6ir + 0.0257' -1.225e~Y' +2.255e-Y" 

where T is in Kelvin and Mi in Joule/mol. 


(2.49-a) 

(2.49-b) 


(2.50-a) 

(2.50-b) 
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2.4.4 Effective Diffusivity of the Pellet; 


The overall mathematical model of the reactor requires the value of the effective 
diffusivity of the catalyst pellet. The effective diffusivity is the combined effect of the 
Bulk diffusivity and the Knudsen diffusivity. These quantities are modeled as follows: 

(a~) Bulk Diffusion : 

The diffusive flux of a component ‘A’ in a mixture of ‘A’ and ‘B’ is represented 
in terms of concentration gradient and according to Pick’s Law is given by: 


(2-51) 

dz 

Here Dab is the bulk molecular diffusion coefficient for component ‘A’ diffusing 
in a mixture of ‘A’ and ‘B’, Ct is the total concentration, Ya is the mole fraction of specie 
‘A’ and dz is the axial distance along the reactor. Considerable attention has been given 
in the literature to predict the mass diffusion coefficient Dab- Assuming ideal gas 
behavior, kinetic theory may be used to show that: 

3 


D 


AB 


P-^T- 


(2.52) 


This relation applies to restricted pressure and temperature ranges only. Bird et. 
al. [1960] Provide detailed discussions of the available theoretical treatments and 
comparisons with experiments. Skelland [1974] and Ried and Sherwood [1966] provide 
more detailed treatments of this subject. Reid et. al. [1987] and Fuller et. al. [1966] also 
discuss the calculation of Dab. A method developed by Fuller et. al. [1966] gives the 
following formula for Dab, 
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l.OlSe'' -r'-'' 


1 


1 


y-s 

^b) ( 2 . 53 ) 

Dab has the units m^/s, P is the pressure in Pa and are diffusional volumes. 
Values of the latter for simple atoms and molecules are given in the table below [Fuller, 
1966]: 


Table 2.1 Atomic and Structural Diffusion Volume Increments 


c 

H 

0 

N 

Cl 

00 

Aromatic Rings 



5.48 

5.69 

19.5 

17.0 

- 20.2 


Table 2.2 Diffusion Volumes of Common Simple Molecules 


H 2 

N 2 

O 2 

Air 

CO 

CO 2 

N 20 

H 2 O 

NH 3 

7.07 

17.90 

16.60 

20.10 

18.90 

26.90 

35.90 

12.70 

14.90 


Note : For systems involving more than two components, the diffusion of each species can 
be treated as if it were diffusing through another single species rather than through a 
mixture by applying a correction factor to Dab [Hill, 1977]. 


(bl Knudsen Diffusion; 

The pores of the catalyst pellet are usually quite small and there are a significant 
number of collisions between the gas molecules and the pore walls. If the pores are 
sufficiently large so that collision between molecules are much more firequent than those 
between molecule and the wall, the mechanism of diffusion is the same as the bulk phase. 
When the mean free path of the diffusing molecules is larger than the mean pore 
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diameter, then the molecules “rattles” down each pore by successive collisions with the 
pore walls. For a single straight cylindrical pore [Hirschfelder et. al. 1954]: 

2 

^K.A - ^ ■ (Pore Radius) • (Mean Thermal Speed of Molecule ' A') (2.54) 


By kinetic theory of gases it can be shown that: 


^K,A ~ ^2 • 


^ j. 




(2.55) 


Where Vgq is the equivalent pore radius in m, T is in °K and the molar mass of the 
diffusion species, Ma is in kg/kmol or g/mol. Frequently the equivalent pore radius is 
taken as the average radius of the pores in the pellet. Alternatively, it may be computed 
from the expression by Satterfield [1963]: 


(BET) p, ’ 


(2.56) 


where BET= Brunauer-Emmet-Teller surface area of the catalyst, (j)p is the 
pellet porosity and is the pellet density. 

(c) Effective Diffusivitv: 

When the Knudsen and bulk diffusivities are significantly different, applicable 
diffusivity is determined by smaller of the two. 

i)=DABorDj,A (smaller of the two) (2-57) 

When the two diffusivities are commensurate, the combined effect is actually 
worse th^ either acting alone. The following equation is adequate for most purposes for 
calculating the applicable diffusivity D , 
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1 


(2.58) 


1 _ 1 

^ ^ AB ^K,A 

The pores in the pellet are not straight and cylindrical; rather, they are a series of 
tortuous, interconnecting paths of varying cross-sectional areas. It would not be feasible 
to describe diffusion within each and every one of the tortuous pathways individually. 
This requires that we modify the diffusion coefficient of the gas in the pores to account 
for: 

(a) The void fraction of the catalyst. 

(b) The tortuosity of the pores. Tortuosity can be defined as the ratio of the distance 
that a molecule travels along the pore(s) and pore interconnections to get between two 
points, to the shortest distance between those points. 

(c) The variation in the cross sectional area of the pores. 

The effective diffiisivity of the catalyst pellet can be measured experimentally. In the 
absence of experimental data, its value may be estimated from the pore size distribution 
and catalyst porosity data given by [Scott Fogler, 1992], 

<j) -cj 

D =1^.D (2.59) 

T 


where, 

is the pellet porosity, cr is the pore constriction factor which accounts for the 
variation in cross-sectional area that is normal to the diffusion and t is pore tortuosity 
factor. 
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2.5 THE HEAT TRANSFER: 


The heat transfer between a fluid flowing through a packed bed and the wall is 
interpreted according to a one-dimensional model characterized by an overall transfer 
coefficient, hfs and by two-dimensional model with one parameter, the coefficient of 
coefficient, hfs and by two-dimensional model with one parameter, the coefficient of 
thermal conductivity kp or two parameters, the effective thermal conductivity Ap and wall 
heat transfer coefficient, ap. 

The mathematical model of the packed bed requires relationships for heat transfer 
processes occurring within the bed. These processes are modeled as described below. 

2.5.1 Thermal Conductivity of Packed Beds: 

The heat transfer Wasch [1971] et. al. studied on heat transfer in packed beds. The 
experimental results for overall heat transfer coefficient hw, effective thermal 
conductivity kp and wall heat transfer coefficient ctw were correlated as a function of 
Reynolds number, flow rates, particle and tube diameter. All correlations of the 
parameters in function of Reynolds number lead to straight lines, which intersect the heat 
transfer coefficient axis. Correlations were also setup between the parameters of the 
various models. 

The energy balance equations for the model require the effective thermal 
conductivity for the packed bed. The equations are, [Dixon and Cresswell, 1979], 
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1 + - 




^wf^T 


1 + 


16 




k 


p J 


(l-<^>) 


\^p J 


and 


(2.60) 


Keff =Kf+Ks 


1 


16 


1 + - 


^ 1 0.1^ 

1 

\^fs^P J 


(1-^ 


\^p J 


(2.61) 


where the subscripts ‘a ’ and ‘r ’ refer to axial and radial respectively, f’ and ‘s ’ refers to 
fluid and solid phase respectively, ‘w ’ refers to the wall and ‘p ’ refers to the single pellet. 
For Rb greater than 50 and Dt / Dp greater than 10, it is reasonable to approximate the 
above equation as, 

(2.62-a) 
(2.62-b) 


K,^=Kf+K 


The radial and axial conductivity of the solid, and has, respectively are given by, 
[Wakao and Kaguei, 1982] 
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= k„^ - 

rs as 




■Ini^ 

Bk, 


(2.63) 


where. 


B=C 


(2.64) 


and C = 1.25 for spheres and 1.4 for crushed particles [Hayes, 1997]. 

The conductivities, V and kaf are correlated using the Peclet numbers for heat transfer, 
given by, [Wakao and Kaguei, 1982] 




0 . 66 ^ 

Re^ Pr 


(2.65) 


(p \ _ GCpDp _ 0.73<f) 0.5 

Jaf ~ , ~ -r. . / 


Re. Pr 


Re,Pr 


( 2 . 66 ) 


where. 


RCj — 


P^pP^s _ PpG 


Bed Reynolds Number 


(2.67) 


Prandtl Number 


( 2 . 68 ) 
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Thus, provided the thermal conductivities of the fluid and the catalyst particles, 
the bed porosity, equivalent particle diameter and the reactor diameter are known, the 
thermal conductivities required for the reactor model can be computed. 

2.5.2 Bed to Wall Heat Transfer Coefficients: 

In the pseudo-homogeneous model, the combined effective heat transfer coefficient at the 
wall, hw.eff, is used. Defining a Nusselt number at the wall. 


Nu 




fw 


where, from [Yagi and Wakao, 1959] 


Nu^ = 0.6 Pr°-'' Re“ 1 < Re^ < 40 


(2.69) 


(2.70-a) 


= 0.2 Pr Re°' 40 < Re, < 2000 
Also, from [Dixon A.G. and Cresswell,1979], we note that, 


(2.70-b) 


^r,ejf ^rf 


(2.71) 


Thus we can find the effective bed to wall convective heat transfer coefficient. 

2.5.3 Overall Heat Transfer Coefficients: 

The heat transfer between the reactor side fluid and the heat transfer fluid 
surrounding the reactor can be divided into three components. 

(a) Convective heat transfer between the reactor fluid and the inside wall, which is 
governed by the inner convective heat transfer coefficient, . 


(b) Conduction through the reactor wall which is proportional to the thermal 
conductivity of the wall, k-^,. 
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(c) Finally, there is heat transfer between the outside reactor wall and the heat 
transfer fluid, governed by the outer heat transfer coefficient, hg. 

The overall heat transfer coefficient, U, combines these three parameters. For a 
reactor tube of outside diameter Do the relation is. 





Dj 1 

D. h. 

in in 


(2.72) 


The outside heat transfer coefficient may be estimated from standard heat transfer 
correlations. The best value for U is still a matter of debate and various correlations have 
been proposed as given in De Wasch and Froment [1972], Froment and Bischoff [1990] 
and Dixon A.G. and Cresswell [1979]. 

2.6 PELLET DEFFUSION EFFECTIVENESS FACTOR: 

In heterogeneous reaction sequence, mass transfer of reactants first takes place 
from the bulk fluid to the external surface of the pellet. The reactants then diffuse from 
the external surface into and through the pores within the pellet washcoat, with the 
reactions taking place only on the catalytic surface of the pores. A schematic 
representation of this two-step diffusion process is shown in Figure 2.2.Thus we observe 
that there are two types of diffusion resistance in the heterogeneous cataljdic reaction: 

(a) External Resistance : This is encountered for diffusion of the reactants or products 
between the bulk fluid and the external surface of the catalyst. The main cause of this 
external resistance is the formation of a concentration boundary layer as the flow takes 
place along the catalyst pellets. This has a tendency to either increase or decrease the 
concentration at the surface relatively to the bulk or the free stream concentration. The 
resulting concentration difference is the driving force for mass transfer. 
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The rate of mass transfer of a component ‘A’ present in the bulk stream can be 
expressed using the mass transfer coefficient [Incropera and DeWitt, 1996], 

Ja=KACka-Cs.a) (2.73) 

Various ways of estimation of the mass transfer coefficient have been proposed. 
Generally, it can be safely calculated from the heat and mass transfer analogy. For heat 
and mass transfer of a single sphere, under forced convection, the following Ranz and 
Marshall [Wakao and Kaguei, 1982] equation is popularly recognized. 



(2.74-a) 

Nu[ = 2 + 0.6 Re° ' Pr“' 


5;i[=2 + 0.6Re°'5c°-"' 

(2.74-b) 

Wm; = C:> Sh[ = 

kf D,, 

(2.75-a) 

Pr= ^ o 5c= ^ 

^AB 

(2.75-b) 


The Nusselt and Sherwood numbers for packed beds are generally higher than 
those predicted from the above equation. This difference reduces as the Reynolds number 
is reduced. For higher Re 5 molds number Wakao and Kaguei [1982] analyzed a large 
amount of packed bed data and proposed the following correlations. 
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iVMi=2 + l.lRe°^ Pr°'-^ 


(2.76-a) 



(2.76-b) 

(2.77-a) 

(2.77-b) 


From the above equations, it is concluded that, in general, to increase the mass 
transfer coefficient and hence the rate of reaction one may either decrease the particle 
size or increase the bulk fluid velocity. 

The relative magnitude of the external mass transfer resistance, as compared to 
the internal pore diffiision resistance may be calculated by the Biot Number for mass 
transfer, defined as. 



Internal mass transfer resistance 
to difihxsion in solid 
external mass transfer resistance 
to convection in a fluid 


(2.78) 


For the present modeling, the external resistance to mass transfer was neglected as 
the Biu value was calculated to be well above 100. 

(h) Internal Resistance ; This is encountered for diffusion of the reactants from the 
e:rtemal pellet surface (pore mouth) to the interior of the pellet and vice versa for the 
^jducts. In a porous catalyst pellet, for the reaction to ocpnr inside the pores, the 
^ipjtants must diffuse into it. As the reactants diffuse into the pellet washcoat, they 
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encounter active sites and some of the reactants react. The reactant concentration 


decreases as the distance into the catalyst pellet increases, giving rise to a concentration 
gradient into the pellet as shown in Figure 2.3. Thus the entire washcoat is not accessible 
to the same reactant concentration. If the diffusion velocity is much lower than the 
intrinsic reaction rate, the concentration gradient can be quite large. When using the 
intrinsic reaction rates it is necessary to know the gas phase concentration at the active 
sites. Therefore it is necessary to have a means of calculating the effective rate of reaction 
in a catalyst with a variable concentration profile. This is achieved by defining 
effectiveness factor, a parameter that accounts for variations in concentration throughout 
the pellet washcoat resulting in the variation in the rate of reaction. 

Average Rate of Reaction for a catalyst pellet 

TJ = 

Rate of Reaction evaluated at surface conditions 


The governing equation for simultaneous diffusion and reaction within the 
catalyst pellet will provide the functional relationship for effectiveness factor. 

Performing the one dimensional, steady state mole balance in spherical co-ordinates for 
isothermal conditions within the pellet and assuming that the diffusive flux based on 
effective diffusivity within the pellet follows Pick’s law, we get. 


[Flow in] — [Flow out] = [Disappearence by reaction] 
or 



(2.79-a) 
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where O is called the Thiele modulus defined as. 



2 1 Dejr 


( 2 . 81 ) 


( 2 . 82 ) 
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All dimension in mm. 


i Qh 

f-A z “4 


L=50.00 H P ^ 

out 



Figure 2.1 : Differential volume element used in developing shell balance for the ID 
packed bed model, non-adaiabtic case. 



Fluid Flow 


Figure2.2: Heat and mass transfer between fluid and solid surface of the pellet. 
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CHAPTER 3 


THE SOLUTION STRETAGY 


The final sets of governing equations for one dimensional, plug flow, non-adiabatic, 
heterogeneous, steady state catalytic reactor have been described in the chapter 2. These 
form a set of coupled ordinary differential equations in fluid and solid phase. The inlet 
data on gas temperature and composition and mass flow rate is known. The transport 
properties of the packed bed along with the stoichiometric and constitutive constraints 
complete the set of equations. The problem although non-linear, is well posed and is 
solved as follows: 

• The basic governing equations in solid phase are solved. 

• Once the solid phase values are known, new values at z=zo+Az in fluid phase are 
solved by Runge-kutta Gill algorithm for solving simultaneous ordinary differential 
equations. 

• The transport properties are handled as a separate subroutine providing the results to 
the main solver. 
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An implicit Thomas algorithm solves the pellet diffusion equation and effectiveness 
factors are then computed in separate subroutines. 




3.1 THE MAIN SOLVER: 

The final set of governing equations may be summarized as follows: 

3.1.1 The Solid Phase Mole Balance Equation: 


Y' =Y^ + 


l.Ox fc ji ^ Pc(^ 






{(3 .5 X X )+ (o .5 x x 17^0 )}x x (l - 




Y^ =Y-^ + 
^co ^co ^ 


{(- 3 -0 X X J+ (1.0 X x )}x Pc x (l " </») 


-CfXk^xa^ 


Y^ = y-^ + 

^ H^O ^ H2O ^ 


- 4-0 X ql" X q x x (l - (/>) 


CfXk^xa^ 


■^COj ^co- 


+ 


-l.Ox rcoX^coX Pc x{l-(j)) 


CfXk^xa^ 


3.1.2 The Fluid Phase Mole balance Equation: 


dz 


|l .0 X X ^ k^xa^xCfX { y^^^ yc^n, )} 

(c^xvj 


< {(3-5xq:^,+0.5xr-*)x<|,-^ xa^xC^x(y^-y^J} 

dz (Cy xvj 


(3.1-a) 


(3.1-b) 


(3.1-c) 


(3.1-d) 


(3.1-e) 


(3.2-a) 


(3.2-b) 


53 





(3.2-c) 


dz 


dYjo ^ {- 4-Q ^ ^c\h, >«i>-k xa^xC^x )} 

dz (C^xvj 


(3.2-d) 


dYgo^ _ { 1 -0 X r^Q X ^ k^xa^xCj-x { y^q^ Y^q^ j| 
dz (C^xvJ 


(3.2-e) 


3.1.3 The Momentum Balance Equation: 



p « d ,- 4 ,’ 


‘^«>^(l-i).(l.75.G) 








(3.3) 


3.1.4 Solid Phase Energy Balance Equation: 


2;=2>- 


fl Pc ^ (^C3//8 ^ ^CO ^ ^^Rco ^^Co)\ 


hj,xa^ 


(3.4) 


3.1.5 Fluid Phase Energy Balance Equation: 


dz v^xpxCp 


(3.5) 
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We know the initial condition of fluid phase i.e. 


J'/(o)=!'/IU : /’{0) = ;>IU : J}(o)=rjU, 


Using these initial conditions, solid phase equations are solved. So now all the 
solid and fluid phase parameters are known before going to next step. 

By applying Runge-Kutta Gill Simultaneous Solver in fluid phase, values at 
z=zo+Az are obtained. Again repeating the above process we obtain these values at 
different sections along z direction. 

3.2 THE PELLET DEFFUSION EQUATION: 

The one dimensional diffusion reaction equation for the catalyst pellet is given by 


d^c, ^ 2 ac, /(c,) q 

dr^ r dr 

@ r = rp=>C, I =Q;@ r = r^^^ 


= 0 


(3.6) 


After incorporating a fictitious cell at the inside boundary (r = and discretizing 
the above equation along with the boundary conditions, using N nodes between rp and 
we get, 


cr*-2c/+c/^' 1 

j — +_ 




r,. 


Cf'-C/"0 



Ar J 


^ J 


= 0 


and 


@ 7 = 1=>C,. = Cs @ j = N 


2Ar 


0 


(3.7) 
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Rearranging equation 3.7, We get, 


bCr +dCl +aCr 



(3.8) 


Where, 


and 


r 


b = 
and 


^Ar 


-K., 


■ 1 2 . _ 

: d = : a = 

Ar' 


1 1 
-+- 


Ar^ r^Ar j 


, _ ‘P 'w 

(iV-1) 




To solve the above set of equations the following steps are followed, 

1. Guess C; at all nodes form j=2 to AT by linear interpolation between the limits 
Cs<C<0. Apply the boundary condition at j=l and 7 =A^+ 7 (fictitious node). 

2. Solve the resulting tridiagonal matrix by Thomas Algorithm. 

3. Update the value of Q and go to step 1 until a converged concentration field is 
obtained. 

4. Find the effectiveness factor of the pellet for CO oxidation using equation 2.80. 


33 COMPARISON OF COMPUTATIONAL STEP SIZE: 

Initially, computational step size of 3 mm was chosen because pellet diameter is 3 
mm and Khandekar S.[2000] used 3mm step size for computations in the homogeneous 
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model. Figures 3. 3-3. 5 show comparison of computational step size of 3 mm, 1mm and 
0.5 mm at r,„=450°C with same mass flow rate on CO and HC conversion efficiency and, 
solid and fluid phase temperatures. 

It can be seen that use of smaller step size results in more precise calculations as it 
updates solid phase parameters at each small step resulting in higher predicted conversion 
rates and temperatures. Comparing values of the parameters on these step sizes, one can 
see that there is a significant difference between calculations for 3mm step size and 1 mm 
step size. But there is very little difference between 1mm and 0.5 mm step size. 
Therefore, further computations are done using 1mm step size as it yields more precise 
results and is computationally more efficient. 

3.4 COMPUTATION FOR CONICAL CATALYTIC CONVERTER: 

Conical catalytic converter is subdivided into five equal parts as shown in figure 
3.6. Inlet velocity is computed from ‘Fluent’ software by giving mass flow rate at a 
particular temperature. Computations are performed in these subsections separately as 
shown in Figure 3.7. For each subsection area is calculated and then inlet velocity for 
subsequent computational steps are calculated conservation of mass. The values of mole 
fractions and gas temperatures are averaged at the outlet of catalyst bed. 
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CO Conversion, 



Axial Distance along the Bed, mm 


Figure 3.3: Effect of computational step size in axial direction along catalyst bed 

on CO conversion, adiabatic heterogeneous model = 450 ° C, mjj= 1 .362 g/s 
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HC Conversion, 



Axial Distance along the Bed, mm 


Bgure 3.4; Effect of computational step size in axial direction along catalyst bed 

on HC conversion, adiabatic heterogeneous model, 7^^^= 450 °C, m j= 1.362 g/s 
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Teiiiperatiire. 



Axial Distance along the Bed, nom 


re 3.5: Effect of computational step size in axial direction along catalyst bed on golid 

and fluid phase temperatures, adiabatic heterogeneous model, T 450 C, mj^= 1.362 g/s 
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All dimensions are in mm. 


Catalyst Bed Thickness, 12mm 



Figure 3.6: Schematic representation of a conical catalytic converter. 



Figure 3.7; Schematic representation of one of the sub-sections of figure 3.6 
along catalyst bed. 





CHAPTER 4 


RESULTS 


The results obtained from the present study are presented in this chapter. The results 

are presented into the following sequence: 

• CO and HC conversion and catalyst bed temperatures with adiabatic heterogeneous 
model, 

• Effect of adiabatic and non-adiabatic conditions on heterogeneous model predictions, 

• Conqjarison of non-adiabatic heterogeneous and homogeneous models, 

• Comparison of heterogeneous and homogeneous models with the experimental results 
on a portable generator, 

• Flow field and pressure variation along flow ft)r a moped the exhaust muffler 
equ^ped with conical catalytic converter using Fluent 5.3,CFD software, 

• Conversion efficiency, temperature and pressure drop for a moped catalytic converter 
with different catalyst bed thickness. 
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• Comparison of performance of conical and cylindrical catalytic converters of equal 
catalyst volume. 

• Comparison of performance of cylindrical converter with two different catalyst 
volume. 

For computation of results input data used are given in Table 4.1. These data were 
obtained experimentally by Khandekar S. [2000] on a two-stroke gasoline engine 
portable generator set. Two different mass flow rates corresponding to lOOW and 600W 

iw 

generator output and the corresponding measured exhaust gas composition as shown in 
Table 4.2 are used. Computations are also done for two different arbitrary inlet gas 
temperatures of 300 and 450°C are used to study response of the model to the inlet gas 
temperature. These temperatures are also typical of genset operation as observed 
experimentally. 

4.1 CO and HC Conversion and Catalyst Bed Temperatures with Adiabatic 
Heterogeneous Model: 

Computations are done for solid and fluid phase conversion of CO and HC and 
outlet temperature with two different mass flow rates corresponding lOOW and 600W 
output. CO conversion is presented on Figs. 4.1 and 4.2, HC on Figs. 4.3 and 4.4. 
Figs. 4.5 and 4.6 show temperatures along the catalyst bed of the converter. It is seen 
fiom Figs. 4. 1-4.4 that for the same inlet temperature, higher conversion of CO and HC 
are obtained at the lower mass flow rates and vise versa. For example, at 300 °C inlet 
temperature, at 1.362 and 1.750 g/s mass flow rates space velocities are 33.53 s"' and 
12.33 s'* respectively. Lower space velocity at 1.362 g/s mass flow rate results in higher 
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residence time of gases in the converter and hence higher rates of reaction and 
conversion. These results are further supported by the higher solid and fluid phase 
temperatures compared to those for 1.750 g/s at mass flow rate of 1.362 g/s (Figs. 4.5 and 
4.6). At inlet temperature of 300 °C, solid phase temperature at the converter outlet 
increased from 429 to 513 °C and fluid phase temperature from 405 to 472°C when mass 
flow rate is decreased from 1.750 to 1.362 g/s. At 450 °C inlet temperature, 
corresponding increase in solid phase and fluid phase temperatures are smaller with 
reduction in only 28 and 1 8 °C respectively. 

Comparing Figs. 4.1 and 4.2 it is seen that for 1.362 g/s mass flow rate when inlet 
temperature is increased from 300 to 450 °C, CO conversion increases from 46.3% to 
92.6%, doubling the conversion. For these computations, inlet composition of the exhaust 
gas was kept the same. HC conversion increases from 21.6 to 100%, which is more than 
four times for the same conditions. HC conversion appears to be more sensitive to 
temperature than CO conversion. For same mass flow rate, catalyst outgas almost 
doubled increases from 513°C to 101 1°C when inlet temperature increases from 300 to 
450 °C(Figs. 4.5 and 4.6). 

Figs. 4.1 and 4.3 further show that CO and HC conversion both for solid and fluid 
phase increase gradually along the catalyst bed at 300 °C inlet temperature. However, at 
higher inlet temperature of 450 °C, CO and HC conversion increase sharply in the first 
half of the catalyst bed and then tapers off in the second half. At high Met temperature, 
reaction rates are significantly higher and hence higher conversion rates and heat release. 
Subsequently, as the concentration of reacting species reduce the conversion rates and 
heat release also reduce in the second half of the catalyst bed. On the other hand, space 
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velocity increases from 33.5 to 42.3s'' when temperature is increased from 300 to 450 °C 
for the same mass flow rate of 1.362g/s. hence, residence time decreases accordingly its 
effect being reduction in conversion and temperature rise. However, the effect of inlet gas 
temperature on rates of reaction being of exponential form far outweighs negative 
influence of reduced residence time which results in very high increase in CO and HC 
conversion as well as solid and fluid temperatures. 

Figs. 4.1 and 4.2 show that at the lower inlet gas temperature the difference 
between solid and fluid phase CO and HC conversion increases along the catalyst bed. 
However, the higher inlet temperature of 450 ^C, this difference first increases slightly up 
to middle of the catalyst bed and then decreases gradually in the second half of the bed. 
The same behavior is seen solid and fluid phase temperatures in Figs. 4.5 and 4.6. This 
also could be attributed to significantly higher rates of reaction in the first half of the 
catalyst bed at 7i„ = 450 °Ccompared to 300 °C on the solid catalyst. The conversion and 
temperature increase at a faster rate at solid surface of the catalyst than in fluid when inlet 
temperature is increased. One may conclude that for higher inlet temperatures, catalytic 
converter of smaller volume could be used. But for practical applications, an engine 
warm-up and inlet temperature at fight loads are also to be catered for and hence has to be 
arrived at a compromise between conversion efficiency and catalyst volume 
t.2 Effect of Adiabatic and Non-Adiabatic Conditions on Heterogeneous Model: 

Figs. 4.7-4.12 show the effect of heat transfer from catalytic converter to the 
urroundings on CO and HC conversion and catalyst bed temperatures predicted by the 
leterogeneous model. Again, the data are shown for two mass flow rates and two inlet 
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gas temperatures as earlier. As expected, CO and HC conversion (Figs. 4.7-4.10) and 
solid and fluid phase temperatures (Figs. 4.1 1 and 4.12) are higher for adiabatic condition 
than non-adiabatic condition. Fig. 4.7 and 4.8 show that heat transfer from the converter 
reduces CO conversion by 1 to 4% in solid phase and 0.8 to 3% in fluid phase. Higher 
differences are obtained at lower mass flow rates and lower inlet temperatures. For 
example, at 1.362 g/s mass flow rate and 300 '’C inlet temperature, reduction in CO 
conversion obtained are 4 and 3% in solid and fluid phase respectively (Fig. 4.7-a). 
Corresponding reduction for 1.750 g/s mass flow rate and 450 °C inlet temperature are 1 
md 0.8% only (Fig. 4.8-b). Similar trends are observed for HC conversion as well. The 
leat transfer from the converter reduces HC conversion by 0 to 5% in solid phase and 0 
0 2.75% in fluid phase. For 1.362 g/s mass flow rate and = 450°C, 100% HC 
onversion is obtained both in solid and fluid phases. 

Similarly, Fig. 4.11 and 4.12 show that heat transfer from the converter reduces 
olid phase temperature by 35 to 18°C and fluid phase temperature by 30 to 15 °C. Higher 
ifferences are obtained at lower mass flow rates and lower inlet temperatures. At 1.362 
mass flow rate and 300 °C inlet temperature, reduction in solid and fluid phase 
anperatures obtained are 35 and 28 °C, respectively (Fig. 4.11 -a). Corresponding 
duction for 1.750 g/s mass flow rate and 450 ^’C inlet temperature are 26 and 23 °C 
'ig. 4.12-b). 

J Comparison of Non-Adiabatic Heterogeneous and Homogeneous Models: 

Predictions of non-adiabatic heterogeneous and homogenous models are 
mpared for two mass flow rates and two inlet temperatures as mentioned above. Figs. 
13 and 4.14 compares CO conversion. Figs. 4.15 and 4.16 HC conversion and Figs. 
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4.17 and 4.18 the fluid phase temperatures. CO and HC conversion, and fluid phase 
temperature predicted by heterogeneous model are higher than the homogenous model. 
Difference in CO conversion is higher for lower mass flow rate (1.362 g/s) than the 
higher mass flow rate (1.750 g/s) for same inlet temperature (Fig. 4.13 and 4.14). For 
example, at inlet temperature equal to 300 ®C, heterogeneous model gives 7.54 and 
9.72% higher CO conversion at 1.750 to 1.362 g/s mass flow respectively. At 7)„ = 
450°C,corresponding difference was 8.62 and 12.4%. 

Similar trends are seen for HC conversion (Figs. 4.15 and 4.16) in which at inlet 
temperature equal to 300 °C, HC conversion for heterogeneous model is 2.00 compared 
o 4.38% by homogeneous model for 1.750 and 4 compared to 10.46% for 1.362 g/s mass 
low rate. At inlet temperature equal of 450 °C and mass flow fate equal to 1.362 g/s, 
.00% HC conversion occurs in heterogeneous model. But heterogeneous model gives 
>2% HC conversion. At mass flow rate of 1.750 g/s conversion in heterogeneous model 
:oes up to 93% compared to 79.4% conversion for homogeneous model. 

Figs. 4.17 and 4.18 compares fluid phase temperature for the two models. At I)„ = 
00°C, heterogeneous model gives higher fluid phase temperature by 66 and 530C for 
.750 to 1.362 g/s mass flow rate at converter outlet. At = 450°C, fluid phase 
anperature at catalyst outlet is higher by 51 and 67°C for 1.750 and 1.362g/s mass flow 
tes respectively. 

.4 Comparison of Heterogeneous and Homogeneous Models with the Experimental 
esults on a Portable Generator: 

The predictions of heterogeneous model are compared with homogeneous model 
id expCTimental results (Khandekar S. [2000]) on a two-stroke 2-kVA gasoline fueled 
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portable genset. Experimental data on CO and HC conversion at different genset output 
and mass flow rates are given in Table 4.2 and gas temperature in Table 4.3. Computed 
data are given in Table 4.4 and are compared with experimental data in Figs. 4.19-4.24. 
CO conversion predicted by heterogeneous model is compared with those of 
homogeneous model and experimental data in Fig. 4.19. Heterogeneous model predicts 5 
to 10% higher conversion (Table 4.4) than the homogeneous model. Up to 400W output 
of the engine, there is a good correlation with experimental data. CO Conversion 
efficiency at outputs higher than 400 W is under-predicted by both the models. 

Fig. 4.20 compares HC conversion efficiency. At lower loads, heterogeneous 
model slightly under-predicts conversion efficiency but at higher loads, it over predicts. 
Homogenous model under predicts throughout the operational range. 

In Figs. 4.21-4.24, catalyst bed temperatures with heterogeneous and homogenous 
models are compared with experimental data at no load, 200W, 400W and 700W outputs, 
•espectively. Heterogeneous model predicts higher bed temperature than homogenous 
node!, although these are still significantly lower than the experimental data. Further the 
ixperimental data show that catalyst bed temperature rise sharply to near maximum at 
list about 17mm downstream of the inlet. Such behavior of heterogeneous model is 
bsorved at 450 °C inlet temperature. However, even at 700W output, the measured inlet 
anperature is only 385 °C and using this temperature model predicts only a gradual 
icrease in gas temperature along the catalyst bed. Experimental exhaust gas outlet 
mperatures at 500W and 700W are 772 and 866 °C respectively. 

It is known that significant oxidation of CO in bulk gas results for temperature 
>ove 700°C and oxidation of HC above 600°C. As the gas-sampling probe is not exactly 
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it the C5onverter outlet but in the mixing volume at the converter outlet, post catal3^ic 
eaction might have occurred showing higher CO conversion and gas temperatures. As 
egards, HC conversion rates of reaction for propane have been used observed in 
lififerential (micro) reactors. The reaction rates for typical gasoline in integral reactors 
ould be different which may account for the difference observed between experimental 
nd predicted values. However, the predicted in CO and HC conversion and gas 
eiqperatures show good correlation with the experimental data. Heterogeneous model 
as shown significant inq)rovement over homogeneous model in correlation with 
)q)erimental observations. 

Flow Field and Pressure Variation for the Exhaust Muffler Equipped with 
!!oiucal Catalytic Converter using Fluent 5.3,CFD Software: 

.5.1 Procedure: 

For application of the model to complicated catalytic converter geometry, inlet 
ichaust mass flow rate, gas condition and gas velocity to different converter sections are 
jquired. Pressure drop across the converter is another important data for ensuring 
HTiinuim adverse effect of the converter fitment in the engine muffler. 

There are many CFD softwares like Fluent, FlowSD, Phoenix etc which offer a 
slatively quicker engineering solution to a variety of flow and heat transfer related 
robtems. These softwares may be effectively used for optimization of design and 
duction of time for prototype development. It is important to note that at this stage, the 
•ea of specific chemical kinetic modeling as required for the catalytic converter system, 
not within the domain of these commercial. The softwares are effective and quick tools 
I get valuable flow information especially in systems involving conq)licated geometry. 
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These softwares numerically solve the governing partial differential equations that 
lescribe the transport processes. The partial differential equations are reduced to 
Jgebraic expressions either by integrating the governing equations (in control volume 
ttethod) or discretizing by ‘Taylor aeries approximations’ (in the finite volume method) 
»ver the discrete sub-domains of the overall physical domain describing the problem, 
lierefijre, a set of grid lines must be specified in order to discretize and create sub- 
dmains. The process of creating this sub-domain is called grid generation. 

In this study grids are generated and boundary conditions are applied using 
Ljambif and for solving the equations, ‘Fluent 5.3’ is used. 

Following steps are followed in ‘Gambit’, 

. Create geometry of the problem concerned. 

. Create the faces so that it can draw grid. 

. Choose the size and type of grid for the geometry. 

, Define the boundary conditions i.e. inlet velocity, outlet velocity, wall, symmetry etc. 

, Define the continuum whether it is solid or fluid. 

Ejqwrt it to ‘Fluent’. 

Following steps are followed in ‘Fluent’, 

Scale the geometry. 

Iiput the material and its physical parameters like density, specific heat, viscosity 
id tiiermal conductivity. 

Choose the model for the flow. 

Define tlK operating condition i.e. value of pressure at any cross-section. 
lE^ut the \^ues of boundary condition. 
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6. If the material is porous like in present case, give the porosity for the domain 
comprising porous material. 

7. Define the velocity vector i.e. direction if the domain of the porous material is not 
perpendicular to flow as in the present case. 

8. Iterate the equation, results are obtained once the convergence criteria is satisfied. 

4.5.2 Geometry of Catalytic Converter and Flow Field (Asymmetric Muffler): 

Fig. 4.26 shows the drawing of the exhaust muffler of a typical moped designed by 
M/s lOCL, R&D Centre, Faridabad. This muffler is fitted with the packed bed catal 5 dic 
converter of the conical configuration with its apex in the down stream direction. The 
muffler has asymmetric exhaust gas inlet. To solve the flow problem on this physical 
domain, rectangular grids of 3 mm size are generated (Fig. 4.27). Fig. 4.28 and 4.29 show 
the velocity contours and velocity vectors respectively. It is seen that the velocity of gas 
through the converter increases as one moves from periphery towards apex. Fig. 4.30 
shows contours of stream function in the muffler. There are also several recirculation 
zones on expansion in the different part of the muffler. 

It may be concluded from this figure that being higher near apex of the cone, the part 
of converter away from cone is relatively underutilized. It is seen that velocity inlet to the 
converter is 2.9 to 2.93m/s in the lower portion of the converter when one moves towards 
^ex of the converter. Inlet velocity in upper portion of the converter is 0.76 to 0.78, 
higher being towards the apex. Fig. 4.31 shows the pressure variation. A pressure 
difference of 1.46 KPa is seen across the muffler. Outlet pressure is atmospheric in all 
cases. 
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Next it is studied whether flow field becomes more uniform if the apex of the 
cone is located towards upstream. Fig. 4.32 shows the schematic of asymmetric exhaust 
muffler with apex of the conical converter in the up-stream direction. Fig. 4.33 is for grid. 
These grids are rectangular with each side is of 3mm. Figs. 4.34 and 4.35 are for velocity 
in contour and vector form respectively. There is higher velocity in the narrow passages 
between the obstruction plates but smaller velocity is in the comers of the inlet portion 
where recirculation of flow occurs. Velocity through the catalyst bed is smaller because 
of obstraction by the pellets. Similarly, Fig. 4.36 and 4.38 are for stream function and 
pressure respectively. Unlike muffler of Fig. 4.30, there is less recirculation of flow in the 
inlet section (4.36), which gives better result. Inlet velocity to the converter is 0.98- 
l.lm/s in the lower portion while in upper portion, it is 0.65-0.67m/s. higher velocity is 
towards apex. There is a pressure difference of 1.63KPa across muffler. 

4.5.2 Geometry of Catalytic Converter and Flow Field (Symmetric Muffler): 

To study the effect of location of exhaust gas inlet on flow field, symmetric inlet 
geometry to the muffler was created as shown in Figs. 4.38 and 4.44. Both versions of 
apex towards downstream as well as upstream were studied. 

Fig. 4.38 gives the drawing of muffler with symmetric inlet and with apex of the 
converter towards downstream direction. Fig. 4.39 shows rectangular grid of 3mm size. 
Figs. 4.40-4.43 show velocity contours, velocity vectors, stream function and pressure 
respectively. Here recirculation of flow in the both comers of the inlet section is observed 
(Figs. 4.40-4.42). It is seen that higher portion of the flow takes place near the converging 
section. This may result in poor conversion as volume of the catalyst near the apex is less 
and also overheating of the pellets in the converging section may take place. So geometry 
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should be modified to avoid this phenomena. Inlet velocity to catalyst bed varies between 
0.78 to O.Slm/s. higher velocity occurs near apex. 2.12 KPa pressure drop occurs across 
the muffler. 

In the same way, Fig. 4.44 shows schematic of exhaust muffler with apex of the 
converter towards upstream direction. Grid is shown in Fig. 4.45 and Figs. 4.46-4.49 are 
for velocity contours, velocity vectors, stream function and pressure drop, respectively. 
As expected, same pattern is seen here. Velocity varies firom 0.77 to 0.80m/s higher being 
near apex of the converter. Here pressure drop of 1.73 KPa is seen which is lower than 
the previous arrangement in which apex of the cone was towards downstream direction. 

For further computations for catalytic conversion, inlet velocity of exhaust gas to the 
converter is taken as uniform as observed in muffler with symmetric inlet velocity 
condition, for example, for mass flow rate of 1.5525g/s, inlet velocity to e converter 
section is taken equal to 0.8m/s for 7^„=300°C. 

4.6 Conversion Efficiency, Temperature and Pressure Drop for Conical Catalytic 
Converter of a Moped with Different Catalyst Bed Thickness: 

The computations are done for these exhaust mass flow rates at different inlet 
temperatures for a fixed exhaust gas composition given in Table 4.6. 

At first, conical catalytic converter of 12mm catalyst bed thickness as designed by 
M/s Indian Oil Corporation Limited, R&D Centre, Faridabad is studied. Five different 
inlet temperature e.g. 300, 350, 400, 450 and 500°C are taken with four exhaust mass 
flow rates typical of a moped at 30, 40, 50 and 65 km/hr as shown in Table 4.5. With 
iiese mass flow rates, inlet velocity to the converter is obtained through Fluent software 
It temperature of 300°C. It is seen that velocity varies slightly in inlet to the converter 
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firom 0.78 to 0.81 at mass flow rate equal to 1.5225 m/s. In the same way for mass flow 
rate equal to 2.349g/s, velocity varies from .98 to l.lm/s. For different mass flow rates 
and inlet temperatures, the gas inlet velocities to the converter section are given in Table 
4.7. 

Figs. 4.50 and 4.51 shows CO and HC conversion at different inlet gas 
temperatures for different mass flow rates 1.5225 to 5.0g/s corresponding moped speed at 
30 to 65Km/hr. At lower mass flow rate, i.e. at lower speed, inlet temperature is also 
lower. So at 300 and 350°C inlet temperature and at 1.5225g/s mass flow rate (30Km/hr 
of vehicle speed), CO conversion is 23.5 and 39% respectively. HC conversion is only 
3.3 and 14.5% respectively for the same condition. However at 450°C inlet temperature, 
CO conversion is 67% and HC conversion is 56%. But in actual operating conditions, at 
lower mass flow rate i.e. at lower speed, inlet temperature Of 450°C may not be obtained. 
It is seen that conversion efficiency especially for HC at practical inlet temperature at 
different mass flow rates is quite low. Thus 12mm bed thickness is not suitable for 
practical purposes. Hence, bed thickness is increased to 15mm and 18mm. For 15mm bed 
thickness, CO and HC conversion, outlet temperature and pressure drop are shown in 
Figs. 4.54-4.57 respectively and for 18mm bed thickness, these are shown in Figs, 4.58- 
4.62. for 15 mm bed thickness, at 30km speed, CO conversion is 33.5 and 55.1% while 
HC is 6.6 and 27.6% at inlet temperatures of 300 and 350°C respectively. Corresponding 
outlet temperatures are 396 and 495°C respectively. Increasing catalyst bed thickness to 
ISmni does not improve HC conversion to the desired label. This shows that there is not 
much improvement in conversion efficiency. So catalyst bed thickness is further 
increased to 18mm. 
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For 1 8mm catalyst bed thickness, CO and HC conversions, outlet temperature and 
pressure drop are shown in Figs. 4.58-4.62 respectively. It is seen that for 1.5225g/s mass 
flow rate and 300 and 350°C inlet temperature, CO conversions are 44 and 68.2% and 
HC conversion is 11.9 and 44% respectively. Although, at higher inlet temperature 
(450°C) and higher vehicle speed (65Rm/hr), CO and HC conversion are significantly 
53.9 and 42% compared to 17.3 and 3.8% at 300*'C inlet temperature. 

It is found from the above discussion that conical converter even with 18nmi 
catalyst bed thickness does not give significant CO and HC conversion especially at 
lower inlet temperature and even at lower vehicle speed. But on the other hand it uses a 
large amount of catalyst volume (192cm^). It should be noted here that catalyst volume 
for 4-stroke engines used are generally equal to swept volume of the engine or lower. For 
the 2-stroke engines, catalyst volume are somewhat higher as higher HC and CO 
emissions are to be converted with adequate catalyst life. This moped is equipped with an 
engine of 75 cm^ swept volume. Very high catalyst volume is not acceptable in practice, 
as the cost will be higher. Conversion efficiency, of equal volume of conical and 
cylindrical shaped converters are compared and discussed in next section. 

4.7 Comparison of Conical and Cylindrical Catalytic Converter Having Equal 
Catalyst Volume: 

Performance of conical converter with 12, 15 and 18mm catalyst bed thickness 
are compared with cylindrical catalytic converter of 117cm^ catalyst volume {Dr — 
50mm, £=60mm) which has the same catalyst volume as the 12mm conical converter at 
different inlet temperatures and exhaust mass flow rates is studied. 
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CO conversions are compared at different inlet temperatures in Figs. 4.63 -4.65 
and HC conversion in Figs. 4.66-4.68. Both CO and HC conversions are significantly 
higher for the cylindrical converter than the conical ones at lower mass flow rate. For 
example, at 1.5225g/s mass flow rate, CO conversion in cylindrical converter is 78.2% 
while conical converters of 12mm, 15mm and 18mm catalyst bed thickness give only 
23.5, 33.2 and 44.1% conversion respectively (Fig. 4.63). At higher inlet exhaust mass 
flow rate (5.0g/s), the difference in conversion between cylindrical and conical converter 
decreases. Cylindrical converter gives 15.5% conversion whereas conical converters give 
10.2, 13.7 and 17.3% for 12mm, 15mm and 18mm catalyst bed thickness respectively. 
Conversion given by conical converter of 18mm catalyst bed thickness is only slightly 
higher than 12mm equivalent cylindrical catalytic converter at mass flow rate of 5.0g/s, 
in spite of catalyst volume of the former being 62%. 

CO conversion at 400°C inlet temperature and different exhaust mass flow rates is 
shown in Fig. 4.65. At lower mass flow rate (1.5225g/s), 100% CO conversion is 
obtained in cylindrical converter while conical converters give only 55.6 to 79.4%, 
highest being for 18mm catalyst bed thickness. At higher mass flow rate (5.0g/s) and 
400°C inlet temperature, CO conversion in cylindrical converter is 35.1% and in conical 
converter, conversion varies from 22.2-40.4% as catalyst bed thickness varies from 12 to 
18mm. CO conversion for Z',>,=350°C, in between conversions for 300 and 400°C inlet 
temperature as expected (Fig. 4.64). 

At 7;„=300°C and mass flow rate of 1.5225g/s, HC conversion is 84.1% for 
cylindrical converter (Fig. 4.66). For conical converters, 3.3, 6.6 and 11.8% HC 
conversions are observed. At higher mass flow rate (5.0g/s), HC conversion for all the 
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converter configurations are very low. For example, at 300°C inlet temperature, 
cylindrical catalytic converter gives 1.6% conversion while conical converter with 
different catalyst bed thickness gives only 0.7, 1.3 and 2.0% conversions. At higher inlet 
temperature (400°C) and lower mass flow rate (1.5225g/s), HC conversion reaches 100% 
in cylindrical converter and 33.8 to 71.4% for conical converter as bed thickness is 
increased. But for the same inlet temperature and at higher mass flow rate (5.0g/s), HC 
conversion in cylindrical converter is 26.6% and in conical converter is 9.17 to 22.1 %. 

Outlet temperatures for cylindrical and conical converters are shown in Figs. 
4.69-4.71. Fig. 4.69 shows that at low mass flow rate (1.5225g/s) and inlet temperature of 
300°C, outlet temperature for cylindrical converter is 823°C and for conical converters 
are only 334to 467°C. But at higher mass flow rate (5.0g/s) outlet temperatures for all 
converter configurations are only 354-362°C. Fig. 4.71 shows outlet temperatures at 
400°C inlet temperature. It is seen that at lower mass flow rate (1.5225g/s), outlet 
temperature for cylindrical converter reaches to 1033®C and for conical converters to 
673and 875°C. At the same inlet temperature but at higher mass flow rate (5.0g/s), outlet 
temperature in cylindrical converter is 604°C and in conical converters varies firom 500 
to 601°C. 

Figs. 4.72-4.73 compares pressure drop of cylindrical and conical convertors. At 
inlet temperature equal to 300*^C and lower mass flow rate (1.5225g/s), cylindrical 
converter gives l.lKiPa pressure drop whereas conical converter gives only 0.4 to O.SKPa 
(Fig. 4.72). At higher mass flow rate (5.0g/s), pressure drop in cylindrical converter 
reaches as high as S.OKPa but in conical converter it is much smaller (of the order of 1.0 
to 1.4KPa). At inlet temperature (400°C) and for the same mass flow rate, pressure drop 
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increases. But a higher increase in pressure drop is observed for cylindrical converter. For 
example, 2.4kPa pressure drop for 1.5225g/s and 7.4kPa for 5.0g/s mass flow rate is seen. 
The corresponding pressure drop for conical converters are only 0.5-0.8kPa and 1.3 to 1.9 
kPa. 

It is seen that cylindrical converter gives much higher converter rates compared to 
conical converter its main disadvantage is that but pressure drop is significantly high. So, 
further study is done on cylindrical converter with lower catalyst volume keeping the 
same diameter (50mm), this is discussed in next section. 

4.8 Comparison of cylindrical converters with two different catalyst volumes: 

Comparison on cylindrical converters with two different catalyst volumes (117 
and 97cm^) is studied. Diameter is kept constant and only length is reduced from 60mm 
to 50mm. Figs. 4.75-4.77 show comparison on CO and HC conversion, outlet 
temperature and pressure drop. At lower mass flow rate (1.5225g/s), and Tin = 400°C, 
both converters give 100% conversion but at 300°C, 98cm^, converter give 60% CO 
conversion compared to 79% for 117cm converter. At 5.0g/s mass flow rate and at 
400°C inlet temperature, corresponding CO conversions are 28 and 37% and at 300°C, 12 
and 16% respectively. HC conversion reaches 100% in both converters at 1.5225g/s and 
reduces by 8% at 5.0g/s mass flow rate. The outlet temperature is lowered by 50°C for 
98cm^ at mass flow rate equal to 5.0g/s. Pressure drop reduces by about IkPa in 98cm^ 
compared to 1 17cm^ cylindrical converter. It is seen that there is still significantly high 
pressure drop even when converter with smaller volume is used. Further geometry 
modifications need to be studied for reducing pressure drop in converter. One such 
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geometry could be a hollow cylinder with relatively small inlet diameter to maintain 
radial fluid flow for lower pressure drop and also smaller catalyst volume. 
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Table 4.1 

List of parameters 


S.No. 

Item 

Dim. 


1. 

Length of the cylindrical catalyst bed . L 

50.50 

mm 

2. 

Internal Diameter of the catalyst bed, Dj 

40.00 

mm 

3. 

External diameter of the catalyst bed, 

42.50 

mm 

4. 

Bulk volume of the catalyst bed 

63.46 


5. 

Weight of the catalyst pellets in the bed 

47.97 


6. 

Number of pellets in random sample of 10 em 

420 


7. 


2015 


8. 1 

Approximate volume of single catalyst pellet 

1.9658e-2 

ESH 

9. 

Net volume of the pellets in the catalyst bed 

39.61 

EiSH 

mm 

Void fraction of the catalyst bed, ^ 

0.375 


11. 

Catalyst: Pt-Rh 

i 

— 

12. 

Catalyst support: Modified cordierite 


— 

13. 

Pellet size 

3.0-3.2 

mm 

14. 

Bulk density 

ro.6 

g/cm^ 

15. 

Pore: Macropores 

~ 

— 

16. 

Pore size 

2 

JL 


Table 4.2 

Experimental data with inlet and outlet values of species and inlet mass flow rate at 

different loads of the genset engine. 


Load 

CO 

HC 

CO 2 

0 

2 

Total 

air + fuel 

Inlet 

Outlet 



Inlet 

Outlet 

Inlet 

Outlet 

% vol. 

PPMCa 

% vol. 

%vol. 

g/s 

Idle 

4.38 

3.82 

5590 



11.18 

4.37 

3.64 

1.128 

100 W 

3.97 

3.28 

4860 

4590 

10.77 

11.68 

3.78 

3.23 

1.362 

200 W 

3.78 

2.47 

4700 

3920 

11.07 

12.74 

3.74 

2.45 

1.367 

300 W 

3.96 

2.15 

5130 



13.06 

3.93 

2.25 

1.458 


3.69 

1.78 

5670 

4050 

1 

10.50 

13.93 

4.30 

1.66 

1.512 


3.36 

1.35 

6120 

4040 

10.30 

14.09 

4.72 

1.80 

1.582 


3.53 

1.30 

6970 

4430 

9.91 

14.16 

5.18 

1.74 

1.750 

700 W 

3.175 

0.74 

7510 

4370 

9.61 

14.40 

5.65 

1.62 

1.943 
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Table 4.3 

Experimental data of temperature along bed at different loads ofthegenset engine 


Load 


Temperature ,"C | 

Tin 

Tn 

T34 

Tout 

Idle 

275 

664 

365 i 

409 

lOOW 

290 

376 

477 

517 

200W 

301 

466 ’ 

499 

525 

300W 

327 

509 

530 

543 

400W 

348 

'601 

599 

657 

500W 

363 

725 

756 

772 

600W 

376 

803 

838 

866 

700W 

382 

824 

846 

866 


where Tn and T34 are temperatures measured from 17mm and 34mm to the inlet of the 
converter. 


Table 4.4 


Computed results of heterogeneous and homogeneous models with experimental 
values using experimental catafytic converter inlet gas data. 


Load 

Tin 

(°C) 

Heterogeneous 

Model 

Homogeneous 

Model 

Experimental 

Values 

CO 

HC 

Tout 

CO 

HC 

Tout 

CO 

HC 

Tout 

(% conversion) 

(‘‘C) 

(% conversion) 

("C) 

(% conversion) 


Idle 

275 

10.10 

1.22 

305 

8.21 

1.01 

297 

12.78 

3.22 

409 

lOOW 

290 

24.22 

3.65 

374 

15.76 

1.81 

338 

17.38 

5.56 

517 

200W 

301 

36.97 

10.22 

442 






525 

300W 

327 

52.41 

33.31 







543 

400W 

348 

54.78 

44.00^ 

650 






657 

500W 

363 

52.65 

46.07 







772 

600W 

376 

48.41 

49.33 







855 

700W 

382 

42.45 

43.10 







866 



















Table 4.5: Measured Exhaust Mass Flow Rate vereus Vehicle Speed 


for a Moped (lOCL Ltd.) 


S. No. 

Vehicle speed 
(km/hr) 

Exhaust mass flow 
rate (g/s) 

1. 

30 

1.5225 

2. 

40 

2.349 

3. 

50 

3.5225 

4. 

65 

5.0 


Table 4.6; Composition of Inlet Gas used for Computation 


CO 

HC 

CO 2 

O 2 

H 2 O 

(% vol.) 

(PPM Cs) 

(%voL) 

(%vol.) 

(%vol.) 

3.96 

5130 

11.07 

1 

! 

3.74 

0.0 


Table 4.7: Calculated Mean Exhaust Gas Velocity to 


Converter Section Using Fluent 


S. No. 

Exhaust 
Mass Flow 
Rate 
(g/s) 

Inlet Gas Velocity 
(m/s) 

Ti„=300^C 




1. 

1.5225 

0.80 

0.87 

0.94 

1.00 

2. 

2.349 

1.00 

1.08 

1.18 

1.22 

3. 

3.5525 

1.20 

1.30 

1.41 

1.51 

4. 

5.000 

1.70 

1.85 

2.00 

2.15 : 
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Axial Distance along the Bed, mm 


Figure 4.1: CO conversion efficiency along catalyst bed, adiabatic heterogeneous model, 
T. =300 9: 

in 



Figure 4.2: CO conversion efficiency along catalyst bed, adiabatic heterogeneous model, 
T. =450 93 

in 
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Figure 4.4: HC conversion efficiency along catalyst bed, adiabatic heterogeneous model, 
Tj^=450 ^ 
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Temperature, 
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CO Conversion, % CO Conversion, % 



(a) 



(b) 

Figure 4.7: Comparison of adiabatic and non-adiabatic heterogeneous model 

(heat transfer to atmosphere from converter) on CO conversion, T jj{-300 
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CO Conversion, % CO Conversion, % 




(b) 

Figure 4.8: Comparison of adiabatic and non-adiabatic heterogeneous model ^ 

(heat transfer to atmosphere from converter) on CO conversion, =450 C 
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HC Conversion, % HC Conversion, % 



15 20 25 30 35 

Axial Distance along the Bed, mm 


Figure 4.9: Comparison of adiabatic and non-adiabatic heterogeneous model 

(heat transfer to atmosphere from converter) on HC conversion, =300 C 
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HC Conversion, % HC Conversion, % 




(b) 

Figure 4.10: Comparison of adiabatic and non-adiabatic o^ 

(heat transfer to atmosphere from converter) on HC conversion, T:, -450 C 
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Temperature, Temperature, 



(a) 



(b) 

?pre 4.11: Comparison of adiabatic and non-adiabatic heterogeneous model (heat transfer 

o 

to atmosphere from converter) on solid and fluid phase temperatures, T = 300 C 
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(a) 



(b) 

Figure 4.12: Comparison of adiabatic and non-adiabatic heterogeneous model (heat transfer 

to atmosphere from converter) on solid and fluid phase temperatures, T j = 450 °C 
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Figure 4.13: Comparison of CO conversion efficiency with non-adiabatic heterogeneous and 
homogneous models, = 300°C 



“1(0=1-362 g/s mjj^=1.750 g/s 

A = Fluid Phase, Heterogeneous Model C = Fluid Phase, Heterogeneous Moldel 

B = Homogeneous Model D = Homogeneous Model 

Figure 4.14 : Comparison of CO conversion efficiency with non-adiabatic heterogeneous and 
homogneous models, = 450°C 
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Figure 4. 15; Comparison of HC conversion efficiency with non-adiabatic heterogeneous and 
horaogneous models, = 300°C 



mjjj=1.362 g/s mjjj=1.750 g/s 

A = Fluid Phase, Heterogeneous Model C = Fluid Phase, Heterogeneous Moldel 
B = Homogeneous Model D = Homogeneous Model 


Figure 4.16; Comparison of HC conversion efficiency with non-adiabatic heterogeneous and 
homogneous models, = 450°C 
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Figure 4.17: Comparison of solid and fluid phase temperatures with non-adiabatic 
heterogeneous and homogneous models, 7^^^= 30CP C 



mjj^=1.362 g/s mjjj=L750 g/s 

A = Fluid Phase, Heterogeneous Model C = Fluid Phase, Heterogeneous Moldel 
B = Homogeneous Model D = Homogeneous Model 

Figure 4.18: Comparison of solid and fluid phase temperatures with non-adiabatic 
heterogeneous and homogneous models, Tjjj= 45CP C 
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Figure 4.19: Comparison of CO conversion (fluid phase heterogeneous model, 
homogeneous model and experimental) 



Figure 4.20: Comparison of HC conversion (fluid phase heterogeneous model, 
homogeneous model and experimental) 
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750 



Figure 4.21: Comparison of predicted and experimental temperature at 
no load with non-adiabatic condition 



Figure 4.22: Comparison of predicted and experimental temperature at 200 W 
engine load with non-adiabatic condition 
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1000 



Figure 4.23: Comparison of predicted and experimental temperature at 400 W 
engine load with non-adiabatic condition 



A = Solid Phase, Heterogeneous Model 
B = Fluid Phase, Heterogeneous Model 
C = Homogeneous Model 
D = Experimenatal 

Figure 4.24: Comparison of predicted and experimental temperature at 700 W 
engine load with non-adiabatic condition 
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Figure. 4.25: Flow chart for solving problems using CFD software. 
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Figure 4.26 : Schematic of asymmetric exhaust muffler with apex of the conical converter 
towards downstream direction 
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Figure 4.27: Grids of exhaust, muffler with apex of the 
converter towards downstream direction 
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Figure 4.28: Velocity contours of exhaust muffler with apex 
of the converter towards downstream direction 


Scale: lmm = 0.5in/s 
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Figure 4.29: Velocity vectors of exhaust muffler with apex 
of the converter towards downstream direction 
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Figure 4.30: Stream function of exhaust muffler with apex 
of the converter towards downstream direction 
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Figure 4.31: Erasure variation of exhaust muffler with apex 
of the converter towards downstream direction 









Figure 4.32: Schematic of the exhaust muffler with apex of the conical converter 
towards upstream direction 
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Figure 4.33: Grids of exhaust muffler with apex of the 
converter towards upstream direction 
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F'igure 4.34: Velocity contours of exhaust muffler with apex 
of the converter towards upstream direction 


Scale: lmm = 0.5iii/s 
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Figure 4.35; Velocity vectors of exhaust muttler with a|>ex 
of the converter towards upstream direction 






Figure 4.36: Stream function of exhaust muffler with apex 
of the converter towards upstream direction 
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Figure 4.37: Pressure variation of exhaust muffler with a|>ex 
of the converter towards upstream direction 
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Figure 4.38: Schematic of symmetric exhaust muffler of with apex of the conical converter 
towards downstream direction 



Grid 

Feb 28, 2001 


FLUENT 5.3 (2d, segregated, ke ) 


Figure 4.39: Grids of exhaust muffler with apex of the 
converter towards downstream direction 
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Figure 4.40; Velocity contours of exhaust muffler with apex 
of the converter towards downstream direction 


Scale; 1mm = O.Sm/s 



Velocity Vectors Cctoed Vekicity Magnitude (m/s) 


Feb 27, 2001 
FLUENT 5.3 (2d, segregated, ke) 


Figure 4.41: Velocity vectors of exhaust muffler with apex 
of the converter towards downstream direction 
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Figure 4.42: Stream function of exhaust muffler with apex 
of the converter towards downstream direction 
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Figure 4.43: Pressure variation of exhaust muffler with apex 
of the converter towards downstream direction 
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Figure 4.44 : Schematic of the symmetric exhaust muffler with apex of the conical converter 
towards upstream direction 
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Figure 4.45: Grids of exhaust muffler with apex of the 
converter towards upstream direction 
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Figure 4.46: Velocity contours of exhaust mufSer with apex 
of the converter towards upstream direction 


Scale: 1mm = 0.5m/s 



Velocity Vectors Colored By Velodty Magnitude (m/s) FLUENT 5.3 (2d, segr^ated,^e) 


Figure 4-47: Velocity vectors of exhsuJBt muffler wth apex 

of tii6 converter towards npstreain direction 




Figure 4.48: Stream function of exhaust muffler with apex 
of the converter towards upstream direction 
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Figure 4.49: Pressure variation of exhaust muffler with apex 
of the converter towards upstream direction 
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-.50: CO conversion efficiency with inlet temperature and exhaust mass flow rate, 
conical converter of 12 mm catalyst bed (catalyst vol.=117cm 



30 40 50 

Vehicle Speed, km/hr 


A:Tj^=300°C D;Tj^ = 450°C 

B:Tjjj=350°C E:Tjj^ = 500°C 

C:T. =400°C 

in 

Figuine 4 . 51 : HC cDeversion efficiency with inlet temperature and exhaust mass flow rate, 
conical converter of 12 mm catalyst bed (catalyst vol.=l 17cm ^ 
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30 40 50 65 

Vehicle Speed, km/hr 


A:T. =300°C D:T =450°C 

in ^ m ^ 

B:T.j^=350°C E:T;^=500°C 

C:t” = 400°C 

in 

Figure 4.52: Converter outlet temperature with inlet temperature and exhaust mass flow rate, 
conical converter of 12 mm catalyst bed (catalyst vol.=117cm^ ) 



A;T. =300°C D:T. =450°C 

IB in _ 

B:T. =350°C E:T. =500°C 
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C:Tj„ = 400°C 

Figiiiie 4*53: Pressure drop across bed with inlet temperature mid exhaust mass flow rate, 
conical converter of 12 mm catalyst bed (catalyst voL=l 17cm^ ) 
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4.54: CO conversion efficiency with inlet temperature and exhaust mass flow rate, 
conical converter of 15 mm catalyst bed (catalyst voL=153cm 
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Figure 4.55: HC conversion efficiency with inlet temperature and exhaust mass flow rate, 
conical converter of 15 mm catalyst bed (catalyst vol.=153cm > 
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6: Converter outlet temperature with inlet temperature and exhaust mass flow rate, 
conical converter of 15 mm catalyst bed (catalyst vol.=153cm^ ) 



A:T. =300°C C;T. =400°C 

in m _ 

B;T.jj=350°C D:Tjjj=450°C 

Figure 4.57: Pressure drop across bed with inlet temperature and exhaust mass flow rate, 
conical converter of 15 mm catalyst bed (catalyst vol.=153cm'* ) 
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Figure 4.58: CO conversion efficiency with inlet temperature and exhaust mass flow rate, 
conical converter of 18 mm catalyst bed (catalyst vol.=192cm 
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Figure 4 59: HC conversion efficiency with inlet temperature and exhaust mass flow rate, 
conical converter of 18 mm catalyst bed (catalyst vol.=192cm ) 
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Figure 4.60: Converter outlet temperature with inlet temperature and exhaust mass flow rate, 
conical converter of 15 mm catalyst bed (catalyst vol.=153cm^ ) 
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Figure 4.6 li Pressure drop across be.d with inlet temperature and exhaust mass flow rate, 
conical converter of 15 mm catalyst bed (catalyst vol.= 153cm ) 
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.62: Comparison of cylindrical (117 cm^) and conical converter of 12mm (117 cm^ ), 
IStom (153 cm^) and 18mm (192 cm^) on CO conversion, Tj = 300 °C 
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Figure 4.63 : Comparison of cylindrical ( 1 17 cm^) and conical converter of 12mm (1 17 cm ), 
15mm (153 cm^) and 18mm (192 cm^) on CO conversion, Tj= 350 °C 
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Vehicle Speed, km/hr 

Figure 4.64: Comparison of cylindrical (117 cm^ ) and conical converter of 12min (117 cm^ ), 
15mm<153 cm^) and 18mm (192 cm^) on CO conversion, Tj = 400 ^C 



3 

Figure 4.65 ; Comparison of cylindrical (1 17 cm^) and conical converter of 12inm (1 17 cm ), 
15mm (153 cm^ ) and 18mm (192 cm^ ) on HC conversion, T. = 300° C 
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Figure 4.68: Comparison of cylindrical (1 17 cm^ ) and conical converter -of 12nmi (1 17 cm^ ), 
15inm (153 cm^ ) and 18mm (192 cm^) on outlet temperature, 300° C 
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Figure 4.69: Comparison of cylindrical (117 cm^) and conical converter of 12mm (117 cm^), 
15mm (153 cm^ ) and 18mm (192 cm^ ) on outlet temperature, Tj = 350° C 
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4.70 : Comparison of cylindrical (117 cn? ) and conical converter of 12nun (117 cm^), 
15inm (153 cm^ ) and 18mm (192 cm^) on outlet temperature, Tjj=400° C 
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Figure4.71: Comparison of cylindrical (117 cm^) and conical converter of 12mm (117 cm^), 
15mm (153 cm^) and 18nun (192 cm^) on pressure drop, T-= 300° C 
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Figure 4.72: Comparison of cylindrical (117 cm^) and conical converter of 12nim(117 cm^), 
15mm (153 cm^) and 18mm (192 cm^) on pressure drop, 350° C 
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Figure 4.73; Comparison of cylindrical (117 cm^) and conical converter of 12mm (117 cm^ ), 
15mm (153 cm^) and 18mm (192 cm^) on pressure drop, T= 400° C 
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A:50inm,50nini(98cm3), T=300°C E : 50nim,50iiun(98 cm^), T =400 °C 

T lu <5 in 

B : 50mm,60inm(l 17 cm^ ) , Tj^=: 300 °C F : 50mm,60mm(l 17 cm^ ) , Tn=400 °C 
C : 50mm,50mm(98 cm^), 350 °C 

D: 50nim,60imn(117 cm^ 350 °C 

Figure 4.74; Comparison of cylindrical catalytic converters with two different catalyst 
volume (98 cm ^ and 1 17 cm ^ on CO conversion 



30 40 50 65 

Vehicle Speed, km/hr 


A:50mm,50mm(98cm^), 300 °C E : 50mm,50mm(98 cm^), T^=400 °C 

B ; 50mm,60nim(l 17 cm^ ) , 300 °C F : 50mm,60mm(l 17 cm^ ) , \=400 °C 

C : 50mm,50mm(98 cm^), 1.^^= 350 °C 

D; 50nun,60mm(117 cm^ ),Tjjj=350 °C 

Figure 4.75 : Comparison of cylindrical catdydc converters with two different catalyst 
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A : 50nmi,50mm(98 cm^), T = 300 °C E : 50imn,50mni(98 cm^), T„=400 °C 

B : 50inm,60mm( 1 17 cm^ ) , 300 °C F ; 50min,60min( 1 17 cm^ ) , 'F„=400 °C 

C : 50mm,50mm(98 cm^), 1^^^= 350 °C 
D : 50imn,60mm( 1 17 cm^ ) , T.^= 350 °C 

Figure 4.76; Comparison of cylindrical catalytic converters with two different catalyst 
volume (98 cm ^ and 1 17 cm ^ on outlet temperature 
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A : 50mm,50mm(98 cmh, \= 300 “C E ; 50mm,50mm(98 cmh, T„=400 °C 

B : 50imn,60mm(l 17 cm ^ ) , \= 300 °C F : 50mm,60inm(l 17 cm^ ) , T;^=400 C 

C ; 50mm,50mm(98 cm^), 350 °C 

D; 50mm,60mm(117 cm^ 350 °C 

Figure 4.77; Comparison of cylindrical catalytic converters with two different catalyst 
volume (98 cm^ and 117 cm^ on pressure drop 




CHAPTER 5 


CONCLUSIONS 


A one dimensional, plug flow, non-adiabatic, heterogeneous, steady state 
mathematical model has been developed for packed bed exhaust catalytic converters for 
SI engine application. Numerical scheme incorporating the simultaneous coupled solution 
of all the five basic governing equations of mole, momentum and energy balance has 
been presented. The model incorporates oxidation of propane, representing commercial 
gasoline and carbon monoxide (CO). Propane (HC) oxidation is taken as two step 
reaction, first to CO, which further oxidizes to carbon dioxide in the second step. 
Although the basic governing equations are developed from the theoretical aspect, the 
supporting expressions for intrinsic rate of reaction of propane and carbon monoxide that 
have been adapted from the literature, renders the overaU model as semi-empirical. The 
predictions of developed model are compared with homogeneous model and 

experimental data. 
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Flow pattern in the muffler of a moped fitted with conical packed bed catalytic 
converter developed by M/s Indian Oil Corporation Limited, R&D Centre, Faridabad is 
studied using a CFD software. The catalytic converter model is applied to this conical 
shaped converter. Predictions of model for the conical catalytic converter are compared 
with cylindrical catalytic converter of equal catalyst volume. 

The following conclusions can be drawn fi'om the present study: 

5.1 Numerical Modeling and Correlation with Experimeatal Data: 

1. The semi-empirical approach for modeling integral catalytic converter is quite 
promising and computationally inexpensive. Due to non-availability of oxidation 
reaction rate data for commercial gasoline, which is a mixture of various hydrocarbons, 
use of propane is taken as the viable option. 

2. Although the formulated problem requires the solution of non-linear coupled 
simultaneous differential equations, the Runge-Kutta algorithm coupled with non-linear, 
one dimensional diffusion equation solved by TDM Thomas algorithm did not have any 
convergence feilure in the solution domain. 

3. Predictions of heterogeneous model on CO conversion correlates well with 
experimental data for a 2-stroke, spark ignition engine generator set. The model, 
however, over predicts HC conversion rates. As regards, HC conversion rates of reaction 
for propane in differential (micro) reactors have been used. The reaction rates for 
typical gasoline in integral reactors could be different which may account for the 
differences observed between experimental and predicted values. 

4. The model predictions on the bed tenqjerature profile for tow inlet gas 
temperature do not correlate well with e 5 q)erimental data. However, at higher exhaust 
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inlet temperature, temperature profile along the bed tends to come closer to the 
experimental trends. Because temperature-sampling probe is not exactly located at the 
converter outlet but in the mixing volume at the converter outlet, post catalytic reaction 
might have occurred showing higher CO conversion and gas temperatures. 

5.2 Study of Flow Field in the Exhaust Muffler: 

Study of flow in the exhaust muffler is done using Fluent 5.1 software. Following 
are the findings of the work. 

It is found that asymmetric exhaust muffler having conical converter with apex 
towards upstream direction has less recirculation of flow though pressure drop is higher 
by 0.17kPa than converter with apex in downstream direction. It is also seen that 
symmetric muffler with both configurations does not give any improvement in flow field 
while pressure drop is almost same as in asymmetric muffler. 

5.3 Performance of Conical Converter; 

1. At first, conical converter with 12mm catalyst bed thickness is studied. But it is 
seen that this converter does not give significant conversion of CO and HC particularly at 
low inlet temperature. 

2. Conical converter of 15mm and 18mm give significant improvement in 
conversion but employs relatively very large catalyst volumes. 

3. Comparison of performance of conical converters of 12, 15 and 18mm bed 
thickness with cylindrical converter having catalyst volume equal to that of 12mm 
conical converter, show that cylindrical converter of aspect ratio of about one gives 
significantly higher conversion than conical converters, but has higher pressure drop in 
the range of 5.0 to 6.0kPa compared to about l.OkPa for converter. 
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5.4 Scope for future work: 

The following points are suggested as a scope of future work: 

1 . Axial conduction and radiation term should be incorporated. Such a model will 
then be able to predict more realistic axial temperature profile. 

2. Ejqperimental investigations should be carried out on integral type reactors to 
augment the available reaction rate data. 

3. It is uncertain whether temperature measured is of pellet surface or of fluid. So a 
more detailed setup is needed for precise temperature measurement. 

4. Further study should be done with different configurations of converter including 
monolith type to minimiz e the pressure drop while obtaining high conversion rates. 
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VI. Appendix A: 


The Run2e - Kutta Algorithm 

We consider the two equations 

with the initial conditions x = and y= y„, when t = /„ . 

^ = <l>{t,x,y) 
at 

Assuming that A t = h, Ax = k, and hy = l, the fourth order Runge-Kutta method gives, 
K=h-f{to,^o^yo\ 

k2 — h • f{tg + "t y A ^ 

1 2 = h ■(p(tg + -j/ZjXjj + Y^p T 

h =h-<p{to +ih,x„+\k2,y,+^l2}, 

K=^-f{K + h>^o+ki,yo+h\ 
l^=h-(j>{t^+h,x^+k2,yo+h), 

x^ =x^+^{k^+ Ikj + 'Ik^ + K) 

6 

3^1 = ^ 0 + 7 (A + 2 ^ 2 +2/3 4-/4) 

6 


132 



In a similar maimer the method can be extended to a system of ‘n ’ simultaneous 
equations. 

Appendix B; 

The Thomas Algorithm 

Step 1: Transformation of the one-dimensional finite difference diffusion equation along 
with the proper boundary conditions to the following matrix form, 


'di 
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0 
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0 ■ 


C, 




'•i 
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Step 2: Computation of all the coefficients di.dz dn , ai,a2, —.an , bi,b2 bn, r i,r2, —r„ 


Step 3: Forward substitution as follows, 


For j = 2 to n, do: 
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d[ihd\J]- 


h\j] 

d\j-\] 


a\j-\] 


Ai]=r[i]- 


d\j-l] 


■[/■-I] 


Thus all the bj are eliminated and we get an upper triangular matrix. 


Step 4: Backward substitution as follows, 


1 

^-d[n] 


For i = (n-1) to 1, step —1, do: 


c[j] = 

4j] 

Thus by solving the above the entire solution Q is obtained. 
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